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SUBJECT:  Instruction  Report  K-81-4,  User's  Guide:  Computer  Program  for 

Design  and  Analysis  of  Cast- In-Place  Tunnel  Linings  (NEWTUN) 


All  Corps  Elements  with  Civil  Works  Responsibilities 


1.  The  subject  user's  guide  documents  a  computer  program  named  NEWTUN  that 
can  be  used  for  designing  and  reviewing  cast-in-place  tunnel  linings.  The 
program  was  originally  developed  by  the  New  England  Division  using  funds 
provided  by  the  Civil  Works  Directorate;  Office,  Chief  of  Engineers.  Although 
this  effort  was  not  directly  part  of  the  Computer-Aided  Structural  Engineering 
(CASE)  project,  this  report  is  being  published  under  the  CASE  tanner  because 
its  subject  matter  is  consistent  with  the  goals  and  objectives  of  the  CASE 
project.  As  is  the  goal  with  all  CASE  projects,  the  intent  is  to  provide 

an  organized,  cost-effective  approach  by  making  available  to  the  structural 
engineer  applicable  computer  programs  ready  for  use  when  the  design  need 
arises . 

2.  Structural  engineers  will  be  readily  able  to  tell  by  the  description 

of  the  programs  and  by  the  examples  given  in  the  report  of  the  applicability 
toward  their  needs.  Detailed  documentation  of  the  programs  may  be  obtained 
from  the  Engineering  Computer  Programs  Library  (ECPL)  of  the  U,  S.  Army 
Engineer  Waterways  Experiment  Station  (WES) ,  Vicksburg,  MS. 

3.  We  strongly  encourage  the  use  of  this  program  where  applicable  throughout 
the  Corps. 
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Preface 


This  report  documents  a  computer  program  called  NEWTUN  that  could 
be  useful  for  the  design  and  analysis  of  cast-in-place  tunnel  linings* 
The  program  is  a  modification  of  an  existing  general-purpose  frame  anal¬ 
ysis  program  called  EFFRAM.  This  report  also  describes  a  procedure  for 
performing  a  tunnel  lining  analysis. 

The  program  was  developed  by  Messrs.  William  J.  Holtham  and  James 
R.  Fay  of  the  U.  S.  Army  Engineer  Division,  New  England,  using  funds 
provided  directly  by  the  Civil  Works  Directorate,  Office,  Chief  of  Engi¬ 
neers  (OCE) .  Mr.  Robert  J.  Smith,  Chief,  Structural  Section,  Civil 
Works  Directorate,  reviewed  the  work  and  was  the  OCE  point  of  contact. 

OCE  also  provided  funds  to  the  Automatic  Data  Processing  (ADP) 
Center,  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES),  to  moni¬ 
tor  the  work  and  publish  this  report.  Dr.  N.  Radhakrishnan,  Special 
Technical  Assistant,  ADP  Center,  and  Mr.  H.  Wayne  Jones,  Computer-Aided 
Design  Group,  ADP  Center,  monitored  the  work  and  reviewed  this  report. 
Mr.  Donald  L.  Neumann  was  Chief  of  the  ADP  Center. 

Although  this  effort  was  not  part  of  the  Computer-Aided  Structural 
Engineering  (CASE)  Project  currently  being  funded  by  OCE,  this  report 
is  published  under  that  banner  because  its  subject  matter  is 
relevant  to  the  CASE  Project. 

Directors  of  WES  during  the  period  of  development  and  the  publi¬ 
cation  of  this  report  were  COL  J.  L.  Cannon,  CE,  and  COL  N.  P.  Conover, 
CE.  Technical  Director  was  Mr.  F.  R.  Brown. 
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Conversion  Factors,  Inch-Fo\md  to  Metric  (SI) 
Units  of  Measurement 


Inch-pound  units  of  measurement  used  in  this  report  can  be  converted  to 
metric  (SI)  units  as  follows: 


Multiply 


BJL 


To  Obtain 


cubic  inches 

feet 

inches 

kips  (1000  lb  force) 

foot-kips  (force) 

kips  (force)  per  cubic  foot 

kips  (force)  per  foot 

kips  (force)  per  square  foot 

pounds  (force)  per  square  inch 

tons  (2000  lb  force)  per  foot 


16.387064 

0.3048 

2.54 

4.448222 

1.355818 

157.087477 

14.593904 

47.880263 

6.894757 

29.187808 


cubic  centimetres 

metres 

centimetres 

kilonewtons 

kilonewton-metres 

kilopascals  per  metre 

kilonewtons  per  metre 

kilopascals 

kilopascals 

kilonewtons 


3 


ELECTRONIC  COMPUTER  PROGRAM  ABSTRACT 


TITLE  OF  PROGRAM 

NEWTUN — Design  and  Analysis  of  Cast-in-Place  Tunnel  Linings 

PROGRAM  NO. 

713  F7  D0063 

PREPARING  AGENCY  | 

U.  S.  Army  Engineer  Division,  New  England  I 

AUTHOR(Sl 

William  J.  Holtham 

James  R.  Fay 

DATE  PROGRAM  COMPLETED 

November  1979 

STATUS  OF  PROGRAM  | 

PHASE 

STAGE 

A.  PURPOSE  OF  PROGRAM 

The  program  computes  joint  deflections  and  member  end  forces  for  continuous 
tunnel  linings  which  are  subjected  to  joint  loads  and  joint  displacements. 
The  tunnel  lining  to  be  analyzed  may  be  founded  in  an  elastic  medium.  The 
program  will  compute  the  external  pressure  on  the  lining  and 
will  perform  iteratioias  to  remove  tension  in  the  analogous  spring  loading, 
if  desired. 


B.  PROGRAM  SPECIFICATJONS 

The  program  is  restricted  to  75  members  and  55  nodal  points. 
The  maximum  nodal  point  difference  is  10. 


C.  METHODS 

Finite  element  techniques  utilizing  the  direct  stiffness  method  of  matrix 
analysis  are  the  basis  for  the  program.  The  elastic  medium  theory  is  based 
on  Winkler's  comcept  of  analogous  springs. 


D.  EQUIPMENT  DETAILS 

The  program  Is  written  in  FORTRAN  IV  and  operates  on  Honeywell  6000  Series 
computers  and  the  Boeing  CDC  CYBER  175  system. 


E,  INPUT-OUTPUT 

Input  allows  computations  for  any  number  of  sections  with  any  number  of  load¬ 
ings.  Output  provides  a  listing  of  descriptive  data  for  the  plane  frame,  joint 
displacements,  external  pressures  with  location  of  compression  and  tension 
zones,  and  member  end  forces.  Input  can  be  in  the  conversational  mode  and  can 
be  manipulated  along  with  output  through  a  terminal  or  through  data  files.  The 
program  has  the  capability  of  generating  nodal  points,  elements,  and  element 

InaAa.  -  -  -  -  - 

F  ADDITIONAL  REMARKS 

1.  Handles  variable  rock  conditions. 

2.  External  springs  can  be  placed  on  any  element. 

3.  External  springs  can  have  one-  or  two-way  action. 

4.  External  springs  can  be  located  on  either  side  of  the  element  by  reversal 
of  the  I  and  J  node  numbers. 

5.  The  program  is  part  of  the  CORPS  system  maintained  by  the  WES  ADP  Center  on 

WES  ,  2205  rtCPLACES  CMC  r  OHM  WHICH  IS  OBSOLETE- 


USER'S  GUIDE:  COMPUTER  PROGRAM  FOR  DESIGN  AND  ANALYSIS 
OF  CAST-IN-PLACE  TUNNEL  LININGS  (NEWTUN) 

Introduction 

1.  The  purpose  of  this  user’s  guide  is  to  demonstrate  the  use  of 
NEWTUN,*  a  computer  program  that  can  be  used  to  analyze  tunnels  with  a 
continuous  cross  section.  NEWTUN  is  a  revision  of  an  existing  general- 
purpose  program  called  EFFRAM**  which  iteratively  analyzes  plane  frame 
structures  on  elastic  foundations.  The  documentation  for  EFFRAM  is  avail¬ 
able  from  the  Engineering  Computer  Program  Library  (ECPL)  of  the  U.  S. 

Army  Engineer  Waterways  Experiment  Station  (WES).  The  documentation  ex¬ 
plains  the  program’s  method  of  solution,  the  stiffness  method  of  struc¬ 
tural  analysis,  and  the  theory  for  beam  elements  on  an  elastic  founda¬ 
tion.  Excerpts  from  the  EFFRAM  documentation  are  included  in  Appendix  A 
to  this  report.  The  main  text  contains  information  pertaining  to  tunnel 
applications  and  data  preparation. 

2.  A  lined  tunnel  is  a  frame  structure  surrounded  by  an  elastic 
medium;  thus,  it  lends  itself  to  analysis  by  EFFRAM.  NEWTUN  has  many  in¬ 
put  options  and  can  handle  any  shape  tunnel  constructed  in  an  elastic  mate¬ 
rial.  Unlike  some  mathematical  models,  NEWTUN  has  springs  over  the  entire 
length  of  an  element  to  simulate  the  subgrade  modulus  of  the  foundation 
material.  (For  tunnels  in  rock,  the  subgrade  modulus  is  referred  to  as 
the  rock  modulus.) 

3.  There  are  a  few  characteristics  of  the  program  that  should  be 
mentioned  for  tunnel  application.  The  first  consideration  is  the  ar¬ 
rangement  of  nodal  points  and  elements.  For  symmetrically  shaped  and 
loaded  tunnels,  it  is  useful  to  model  only  half  the  structure  using 
appropriate  end  restraints.  The  example  problems  in  Appendix  B  demon¬ 
strate  this  technique.  Of  course,  the  entire  structure  can  be  modeled; 

*  NEWTUN  is  designated  X0055  in  the  Conversationally  Oriented  Real-Time 
P^rogram-Generating  System  (CORPS)  library.  Three  sheets  entitled 
’’PROGRAM  INFORMATION”  have  been  hand-inserted  inside  the  front  cover 
of  this  report.  They  present  general  information  on  the  program  and 
describe  how  it  can  be  accessed.  If  procedures  used  to  access  this 
and  other  CORPS  library  programs  should  change,  recipients  of  this  re¬ 
port  will  be  furnished  a  revised  version  of  the  ’’PROGRAM  INFORMATION.” 

**  EFFRAM  is  program  X0022  in  the  CORPS  library. 


however,  it  will  be  necessary  to  have  one  node  restrained  in  either  the 
global  X  or  Y  direction  or  else  erroneous  displacements  will  occur. 

Using  more  nodes  and  elements  will  increase  the  accuracy  of  the  results. 
The  program  has  node  generating  capability;  however,  at  present  circular 
segments  cannot  be  generated.  Certain  load  types  should  be  used  for 
various  loading  conditions.  For  hydrostatic  loads  on  curved  linings, 
nodal  point  loads  should  be  used  by  replacing  the  resultant  load  with 
its  global  X  and  Y  components.  Element  loads  and  loads  that  vary 
linearly  or  uniformly  can  be  used  on  straight-line  segments.  The  use  of 
llneraly  distributed  element  loads  on  curved  sections  should  be  avoided. 

4.  NEWTUN  does  have  some  limitations.  Any  discontinuous  lining 
(e.g.,  a  precast  segmental  concrete  lining)  cannot  be  modeled  because 
it  is  not  possible  to  insert  a  hinge  within  a  section.  It  sho»ild  be 
noted  that,  when  a  precast  tunnel  lining  is  grouted  after  being  placed, 
it  can  act  as  a  continuous  lining  until  it  is  subjected  to  loads  that 
may  crack  the  grout  in  a  segment  joint. 

General  Description  for  Tunnel  Analysis 

5.  This  section  will  highlight  the  information  in  Engineer  Manual 
1110-2-2901  (Headquarters,  Department  of  the  Army  1978)  that  is  essen¬ 
tial  for  modeling  continuous  tunnel  lining.  Discussion  will  be  confined 
to  tunnels  in  rock.  Two  variables  that  have  a  profound  effect  on  the 
usefulness  of  this  program  for  tunnels  in  rock  are  the  rock  load  on  the 
lihing  and  the  rock  modulus.  These  items  are  discussed  in  EM  1110-2- 
2901.  Several  factors  to  be  considered  are  presented  in  Chapter  3, 
"Design,”  of  the  EM.  It  is  suggested  that  the  user  become  familiar  with 
the  contents  of  this  reference  and  he  exercise  keen  engineering  judgment 
when  analyzing  a  tunnel  lining. 

6.  There  are  many  factors  to  consider  when  studying  a  tunnel 
problem.  A  thorough  examination  of  the  geologic  conditions  is  important 
if  a  representative  model  is  to  be  made.  Rock  conditions  are  described 
in  EM  1110-2-2901  in  paragraph  2-4,  "Interpretation  of  Geologic  Data." 
Rock  loads  based  on  the  general  rock  condition  are  outlined  in  Table 
3-1,  "Estimate  of  Rock  Load,"  of  the  EM.  These  values  should  be  used 
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onx/  for  preliminary  designs;  the  final  design  should  incorporate  values 
obtained  from  District  geologists.  Unlike  many  other  types  of  struc¬ 
tures,  there  are  no  standard  loading  conditions  to  be  applied  to  tunnels 
An  explanation  of  the  possible  loads  acting  on  a  lining  is  found  in 

EM  1110-2-2901  in  paragraph  3-7,  "Permanent  Tunnel  Linings,"  The  list 

of  general  design  conditions  in  subparagraph  3-7. b(8),  "Design  Guidance, 
is  a  guide  and  it  is  likely  that  more  load  conditions  may  be  induced  for 
a  particular  tunnel. 

7.  NEWTUN  is  an  attractive  tunnel  analysis  program  because  it  has 

an  option  that  permits  simulating  the  actual  rock  mass  conditions  by  in¬ 

putting  the  rock  modulus.  The  rock  modulus  used  in  NEWTUN  is  applied 
over  an  entire  element,  not  just  at  the  nodes.  (See  Appendix  A  for  the 
theory  of  the  program  and  the  elastic  supports.)  This  feature  results 
in  an  accurate  simulation  of  the  in  situ  conditions  provided  an  adequate 
rock  modulus  is  determined.  An  extensive  treatment  of  selecting  rock 
modulus  is  presented  in  EM  1110-2-2901  in  subparagraph  3-7b(6),  "Method 
for  Determining  Value  of  Foundation,"  and  in  paragraph  7.2.2,  "Supple¬ 
mentary  Comments  on  Significant  Geological  and  Mechanical  Properties." 

It  is  possible  under  some  circumstances  to  have  a  rock  modulus  that 
varies  around  the  lining.  There  are  other  computer  programs  and  theo¬ 
ries  to  assist  the  user  in  such  situations. 


Preparation  of  Data 

8.  The  frame  geometry  is  defined  by  numbering  each  nodal  point 
and  each  frame  member.  No  restrictions  are  imposed  on  the  pattern  of 
nodal  point  or  member  numbering  systems.  However,  the  total  number  of 
nodal  points  cannot  exceed  55,  with  the  maximum  difference  between  mem¬ 
ber  nodal  points  of  10.  The  total  number  of  members  is  restricted  to 
75. 

9.  All  nodal  point  coordinates  are  given  in  the  global  system 
(Figure  1)  and  must  be  positive.  Assuming  the  paper  acts  as  the  global 
plane,  the  X  axis  is  horizontal  and  the  Y  axis  is  vertical.  The  input 
sign  convention  for  loads  calls  for  positive  horizontal  forces  to  the 


V 


Y 

i 


+  MOMENT 

0 


+  FORCE 
IN  Y 

i 


+  FORCE 
IN  X 


( - -  X 

Figure  1.  Global  coordinate  system 

right,  positive  vertic*”*  forces  upwards,  and  positive  moments 
counterclockwise . 

10.  The  boundary  conditions  are  then  determined  in  global  coordi¬ 
nates.  The  type  of  restraint  provided  at  each  nodal  point  for  the  X 
direction,  Y  direction,  and  rotation  is  described  using  a  three-digit 
code  (see  Figure  2).  A  zero  (0)  indicates  that  that  node  is  free  to 
move  in  the  specified  direction  and  that  a  load  is  to  be  specified  at 


Y 


Figure  2.  Typical  nodal  point  restraint  codes 
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the  node*  A  one  (1)  indicates  that  the  node  is  fixed  in  the  correspond¬ 
ing  direction  and  that  a  displacement  and/or  rotation  is  to  be  specified 
at  the  node, 

11.  The  program  allows  both  nodal  point  loads  and  element  loads 
to  be  entered  as  data.  Nodal  point  loads  include  a  force  in  the  X  and 
Y  directions  and  a  moment  about  the  node.  The  sign  convention  is  shown 
in  Figure  1.  The  element  loads  are  of  the  following  types: 


a-  Type  1 
b.  Type  2 
£.  Type  3 
Type  4 

as  shown  in  Figure  3. 


concentrated  load, 
uniformly  distributed  load, 
concentrated  moment, 
triangular  distributed  load. 


TYPE  1 


TYPE  4 


TYPE  2 


TYPE  3 


ELEMENT 


program*  or  by  being  typed  into  a 

data 

file  which  is  read  by  the  program 

If  a  data  file  is  used,  it  should 

have 

line  numbers  at  the  time  it  is 

read  by  the  program. 

Section  title  1 

A. 

Appears  with  section  properties 

Section  title  2 

A. 

Appears  with  section  properties 

Section  title  3 

A. 

Appears  with  section  properties 

Specification  line  1 

A. 

Number  of  nodal  points 

B. 

Number  of  elements 

C. 

Number  of  loading  conditions 

Nodal  point  data 

A. 

Nodal  point  number 

(see  paragraphs 

B. 

Restraint  code** 

13-15  for  genera¬ 

C, 

X  coordinate 

tion  of  nodal 

D. 

Y  coordinate  ^ 

point  data) 

E. 

Displacement  in  X  direction^ 

F. 

Displacement  in  Y  direction  ^ 

G. 

Rotation  about  node  (radians) 

Element  data  (see 

A. 

Member  number 

paragraphs  16-18 

B. 

Node  number  on  end 

for  generation  of 

C. 

Node  number  on  **J**  end 

element  property 

D. 

Smallest  cross-sectional  area  of 

data) 

member  2 

E. 

Young’s  modulus  (force /length  ) 

F. 

Foundation  spring  constant 
(force/  length^) 

G. 

Moment  of  inertia  of  cross 
section 

H. 

Relative  flexural  stiffness  of  . 
end  I  when  end  J  is  fully  fixed*^ 

*  See  Appendix  B  for  example  problems  demonstrating  the  conversational 
mode  of  data  input. 

**  This  is  a  three-digit  code.  The  first  digit  is  for  the  X  direction, 
the  second  digit  is  for  the  Y  direction,  and  the  third  is  for  rotation. 

=  0  if  the  load  is  to  be  specified  or  displacement  is  unknown 

=  1  if  restraint  is  to  be  specified 

t  This  is  specified  only  if  the  restraint  code  equals  one  in  that 
direction. 

3 

tt  Foundation  spring  ■  modulus  of  subgrade  reaction  (force/length  )  x 
width  of  element  (length). 

/  For  a  prismatic  member,  enter  0  and  the  computer  will  compute  the 
value • 


Element  data 

I. 

Relative  flexural  stiffness  of 

(continued) 

end  J  when  end  I  is  fully  fixed* 

J. 

Flexural  carry-over  factor  from 
end  I  to  end  J* 

Loading  condition 

A. 

Appears  with  loading  and  compu¬ 

title  line 

tation  output  only 

Specification  line  2 

A. 

Number  of  nodal  point  load  en¬ 
tries  that  follow 

B. 

NITER:  =  0  ,  iterations  to  re¬ 
move  tension;  =  -1  ,  no  itera¬ 
tions  to  remove  tension 

Nodal  point  load  data 

A. 

Nodal  point  number 

B. 

Force  in  X  direction  (global)** 

C. 

Force  in  Y  direction  (global)** 

D. 

Moment  about  node  (global)** 

Element  load  data 

A. 

Element  number 

(see  paragraph  19 

B. 

Load  type  code; 

for  generation  of 
element  loads) 

=  1  concentrated  load 
=  2  uniformly  distributed  load 
=  3  concentrated  moment 
-  4  triangular  distributed  load 

C. 

Magnitude  of  load: 

for  load  types  1  and  3,  total 
force  on  element 
for  load  type  2,  force  per  unit 
length  on  element 
for  load  type  4,  force  per  unit 
length  on  element  at  weighted 
end 

D.  Location  code: 

for  load  type  1,  =  global  dis¬ 
tance  from  I  node  to  forcet 
for  load  type  2,  =  0 
for  load  type  3,  =  distance 
along  element  from  I  node  to 
moment 

_ (Continued) _ _ _ _ 

For  a  prismatic  member,  enter  0  and  the  computer  will  compute  the 
value. 

This  is  specified  only  if  the  restraint  code  at  the  node  equals  ^ero 
in  that  direction. 

Use  global  X  distance  with  loads  in  global  Y  direction  and  global  Y 
distance  with  loads  in  global  X  direction. 


D.  Location  code  (continued): 

for  load  type  4, 

=  0  if  node  I  is  weighted 
=  1  if  node  J  is  weighted 

E.  Direction  code: 

for  load  types  1,  2,  and  4 
force  direction  code 

=  1  for  Y  axis 
=  2  for  X  axis 

for  load  type  3=0 

Specification  line  3  A.  Indicates  end  of  member  input 

(At  this  point,  repeat  the  procedure  from  the  loading 
condition  title  line  if  more  than  one  loading  condition 
is  to  be  analyzed;  or  repeat  the  procedure  from  section 
title  1  if  more  than  one  section  is  to  be  analyzed.) 

End  of  data  A.  Type 

$END 

to  terminate  input  data 


Element  load  data 
(continued) 


Nodal  point  data  generation 

13.  Restrictions .  The  nodal  point  data  lines  must  be  in  in¬ 
creasing  order  of  nodal  point  numbers  (1  to  n) . 

14.  Generation.  Nodal  point  data  generation  occurs  when 

the  difference  between  the  nodal  point  number  of  the  variable  read  and 
that  of  the  previous  variable  is  greater  than  one.  It  should  be  noted 
that: 


_a.  The  restraint  code  and  specified  displacement  of  any 
generated  nodal  point  is  always  equal  to  zero. 

b.  The  values  of  the  X  or  Y  coordinates  of  the  generated 
node  may  increase  or  decrease.  The  increment  is  com¬ 
puted  by  taking  the  difference  in  the  coordinates  given 
and  dividing  by  the  number  of  elements  to  be  generated. 

Element  property  data  generation 

15.  Restrictions; 

_a.  The  element  data  lines  must  be  in  order  of  increasing 
element  number  (1  to  n) . 

If  the  subgrade  modulus  is  to  be  generated  in  a  linear 
distribution,  the  length  of  each  element  must  be  the 
same. 
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* ' 


16.  Generation.  Generation  occurs  when  the  difference  in  succes¬ 
sive  element  numbers  is  greater  than  one.  It  should  be  noted  that: 

ai.  The  nodal  point  number  of  the  I  node  and  J  node  may  in¬ 
crease  or  decrease  between  the  numbers  given.  The  in¬ 
crement  is  computed  by  taking  the  difference  between 
the  numbers  given  for  the  I  and  J  nodes. 

The  area,  modulus  of  elasticity,  moment  of  inertia,  stiff¬ 
ness  coefficients,  nnd  carry-over  factors  are  set  to  those 
of  the  previous  element . 

The  modulus  of  subgrade  reaction  is  generated  linearly 
from  the  previous  value  to  the  given  value.  (Refer  to  the 
restriction  in  subparagraph  15b  above.)  The  generated 
value  is  an  average  for  the  element .  If  the  modulus 
distribution  is  triangular,  the  value  given  should  be  the 
largest  value  along  the  element.  That  value  will  be  the 
average  for  the  element. 

Element  load  data  generation 
(load  types  2  and  4) _ 

17.  Restrictions: 

The  element  number  must  be  greater  than  100.  The  number 
represents  the  beginning  element  and  the  final  element 
the  load  is  to  be  generated  over. 

jb.  The  elements  need  not  be  the  same  length  since  the  lengths 
are  computed  from  the  coordinate  values. 

c*  For  load  type  4,  the  first  part  of  the  number  represents 
the  element  from  which  the  coordinate  of  the  I  node  is 
taken.  The  last  part  represents  the  element  from  which 
the  coordinate  of  the  J  node  is  taken. 

_d.  For  load  type  4,  the  I  and  J  nodes  must  be  on  the  same 
ends  for  all  elements. 


Program  Diagnostics 

18.  The  following  are  some  programmed  stops  inside  the  program 
with  appropriate  messages  to  explain  the  errors  which  are  encountered: 

***  DATA  ERROR  -  LTYPE  INCORRECT,  LOAD  OMITTED  ELEMENT  _  *** 

The  load  type  code  entered  is  a  number  other  than  1,  2,  3,  or  4.  This 
load  is  omitted  from  any  computations. 
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ERROR  -  THE  NODAL  POINT  DIFFERENCE  AT  BAR  #  IS  = 

MAX  DIFF  =  10 

The  nodal  point  difference,  which  determines  the  band  width  of  the  total 
structure  stiffness  matrix,  is  greater  than  10. 

***  ITERATIONS  HAVE  NOT  CONVERGED,  MAX  =  10  *** 

The  program  is  limited  to  10  iterations  to  remove  tension  in  the  elastic 
springs,  and  is  dependent  on  the  variation  of  the  limits  of  integration. 
If  the  limits  vary  by  more  than  0.05  length  units,  an  iteration  is 
required. 
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Appendix  A;  Excerpts  from  the  Documentation  of 
Computer  Program  EFFRAM* 

Introduction 

1.  The  "Plane  Frame  Analysis  by  Direct  Stiffness"  computer  pro- 

'  aru  acquired  from  the  Rock  Island  District  (Ashton  and  Meyers  1977)  and 
subsequently  revised  by  the  New  England  Division  provided  a  production- 
oriented  program  for  matrix  analysis  of  frar^  ^.iructures.  However,  anal¬ 
ysis  is  limited  to  those  structures  restrained  in  space  by  various  types 
of  supports  (i.e.,  fixed,  pin,  roller,  evc.)*  Structural  problems  are 
often  encountered  where  the  structure  is  founded  on  an  elastic  material, 
The  foundation  reaction  is  dependent  h  the  stiffness  of  the  struc¬ 

ture  and  the  elastic  foundation.  Eduard  Winkier  first  presented  a  theo¬ 
retical  approach  to  the  problem  of  elai^tic  foundations  in  1867.  The 
idealized  foundation  is  analogous  to  an  infinite  number  of  independent 
springs,  where  the  elastic  constant  corresponds  to  the  modulus  of  sub¬ 
grade  reaction  of  the  foundation.  The  reaction  stress  is  proportional 
to  the  deflection.  The  "Beams  on  Elastic  Foundations"  program  written 
by  Professor  John  T.  Christian  of  MIT  applied  the  elastic  theory  to  pris¬ 
matic  beam  elements. 

2.  EFFRAM  combines  and  expands  the  analysis  capabilities  of  each 
program  into  a  single  program  that  can  analyze  frame  structures  on  an 
elastic  foundation.  In  this  program,  no  element  need  be  on  an  elastic 
foundation,  but  the  option  has  been  provided. 

General  Description 

3.  The  plane  frame  is  composed  of  members  lying  in  one  plane  and 
having  axes  of  s)nnmetry  in  that  plane.  The  joints  between  the  members 
are  considered  as  rigid  connections.  Forces  and  displacements  acting  on 


*  Complete  documentation  of  EFFRAM  is  given  in  Holtham  and  Fay  and  is 
available  from  the  Engineering  Computer  Program  Library  (ECPL)  of  the 
U.  S.  Army  Engineer  Waterways  Experiment  Station. 


the  frame  are  in  the  plane  of  the  frame.  Each  nodal  point  (joint)  be¬ 
tween  elements  (members)  has  three  potential  displacements.  These  three 
degrees  of  freedom  consist  of  displacements  along  either  of  the  global 
axes  or  rotation  about  the  joint.  For  any  one  element,  there  are  two 
nodal  points,  each  having  three  degrees  of  freedom.  These  displacements 
are  related  to  the  applied  member  forces  by  the  stiffness  of  the  element 
and  by  the  stiffness  of  the  elastic  foundation. 

4.  A  frame  member  stiffness  matrix  (six  simultaneous  equations, 
each  consisting  of  a  force-stif fness-displacement  relationship)  can  be 
established  for  each  member  in  the  structural  system.  The  superposition 
of  the  individual  member  stiffness  matrix  provides  the  overall  or  total 
structure  stiffness  matrix.  The  solution  follows  once  the  total 
stiffness  matrix  is  established. 

Method  of  Solution 

5.  The  method  of  solution  used  in  EFFRAM  begins  with  computation 
of  the  properties  of  tne  individual  member  stiffness  matrix  in  the  local 
coordinate  system  of  the  member.  If  the  member  is  founded  on  an  elastic 
foundation,  the  foundation  stiffness  matrix  is  computed  and  added  to  the 
member  stiffness  matrix.  The  new  matrix  is  transferred  to  the  assumed 
global  coordinate  system  using  a  transformation  matrix.  The  member 
stiffness  in  global  coordinates  must  now  be  modified  to  account  for 
specified  boundary  conditions  consisting  of  specified  displacements. 

Since  finite  element  techniques  allow  the  addition  of  the  stiffness 
matrices,  the  displacement  restraint  need  not  be  specified  in  the  direc¬ 
tion  of  the  elastic  foundation  (perpendicular  to  the  element).  The 
modified  frame  element  stiffness  matrices  are  next  added  to  form  the 
total  frame  structure  stiffness  matrix  which  relates  the  structure  dis¬ 
placements  and  the  applied  loads.  This  matrix  is  banded  and  symmetrical. 
It  is  necessary  to  store  and  operate  on  the  band  using  a  modified  ver¬ 
sion  of  Gaussian  elimination  to  determine  the  displacements.  From  this, 
the  member  forces  and  reaction  stress  of  the  elastic  foundation  are 
determined. 


6.  In  the  analogous  springs  of  the  Winkler  model,  the  foundation 
properties  are  considered  identical  in  tension  and  compression,  thus 
implying  that  the  foundation  can  support  tensile  stress.  For  most  cases 
this  is  not  a  real  assumption.  In  order  to  take  into  account  the  effects 
of  beam  uplift,  an  iterative  process  was  developed  whereby  only  those 
areas  of  the  beam  element  in  compression  are  used  in  computing  the 
foundation  stiffness  matrix. 

7.  The  foundation  stiffness  matrix  is  composed  of  a  set  of  inte¬ 
gral  equations  in  terms  of  a  dimensionless  parameter.  To  find  each  term 
of  the  foundation  matrix,  the  equations  are  evaluated  over  the  limits  of 
the  compression  areas  in  terms  of  this  parameter.  In  the  first  analysis, 
the  springs  are  on  the  total  element.  If  tension  areas  occur,  iterations 
are  performed  until  the  limits  no  longer  change  (+  0.05  units). 

8.  Unusual  circumstances  may  arise  where  the  foundation  will  have 
two-way  action  (tension  or  compression).  The  program  has  an  option 
which  will  allow  this  action. 


Appendix  B;  Example  Problems 

1.  Two  example  problems  are  presented  in  this  appendix  to  illus¬ 
trate  the  use  of  NEWTUN.  Both  examples  model  the  tunnel  of  the  Park 
River  Auxiliary  Conduit,  Hartford,  Conn.  The  tunnel  is  circular  with  a 
22-tt  inside  diameter  and  a  concrete  liner  13  in.  thick.  The  modulus  of 
elasticity  of  the  liner  is  519,119  ksf,  and  the  surrounding  rock  has  a 
rock  modulus  of  1000  kef.  A  1-ft  length  of  tunnel  is  used  for  the 
analysis . 

2.  Comparisons  were  made  between  NEWTUN  and  the  method  outlined 
in  Engineer  Manual  1110-2-2901  (Headquarters,  Department  of  the  Army 
1978)  and  results  from  another  program  called  TUNNEL  (Orenstein  1973). 
Both  comparisons  showed  favorable  results. 

Example  Problem  1 


3.  A  vertical  rock  load  of  1  kip/ ft  is  the  only  loading  case  used 
in  this  example  (Figure  Bl) .  Included  in  this  example  are  the  loading 
ROCK  LOAD  =  1  kip/ft 


CONCRETE  LINING: 
r  =  4,000  psi 

E  =  519,119  kef 

2 


diagram,  typical  computer  input  and  output,  and  results  from  NEWTUN. 
Figure  B2  shows  the  computer  model;  Figure  B3  shows  the  rock  modulus 
computation.  Table  B1  gives  a  sample  execution  of  the  program  with 
conversational  input.  Table  B2  shows  the  output  from  this  run. 

4.  The  vertical  1-kip/ft  rock  load  is  modeled  as  a  set  of  concen¬ 
trated  loads  applied  at  the  nodal  points.  This  is  the  normal  method  of 
converting  a  distributed  load  to  a  set  of  concentrated  (lumped)  loads; 
i.e. , 


The  applied  load  at  any  node  is  the  sum  of  the  load  half  way  to  each 
adjacent  node.  For  example,  for  node  18, 

left  side  =  1/2  [l  kip/ft  (2.00  -  0.0)]  =  1.00  kip 
+  right  side  =  1/2  [l  kip/ft  x  (3.95  -  2.0)]  ^  0.97  kip 

TOTAL  1.97  kip 


For  node  19, 


left  side  =  (no  load  to  left)  =  0.0  kip 

+  right  side  =  1/2  [1  kip/ft  x  (2.00  -  0.0)]  =  1.00  kip 

TOTAL  1.00  kip 

The  remaining  nodal  loads  can  be  calculated  using  this  same  procedure. 
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GIVEN: 

Elastic  Modulus  =  E  =  1000  ksi  ]  obtained 

^  from  lab 

Rock  Quality  Designation  =  RQD  =  0.85  j  tests 


FIND: 

Rock  Modulus  =  k 

from  Figure  7.2,  page  K-123  of  EM  1110-2-2901, 

the  Modulus  Reduction  Factor  -^  =  0.22 

^r 

where  E  =  elastic  modulus  of  rock  mass 

E  =  elastic  modulus  of  rock  material 
r 

therefore,  if 

^  =  0.22 
^r 

E  =  0.22E 

r 

=  0.22(1000  ksi) 

=  220  ksi 

use  Equation  3-21,  page  3-49  of  EM  1110-2-2901, 
to  obtain  k 


k  = 


(1  +  v)d 


where 


V  =  Poisson's  ratio  =  0.3 
d  =  borehole  diameter  =  24.17 


220 


^  ■  (1  +  0.3)(24.17  X  12) 


h  7?) 


=  1008  kef 
use  k  =  1000  kef 


Figure  B3.  Rock  modulus  computations  for  example  1 
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Table  B1 

Sample  Execution  of  NEWTUN  with  Conversational  Input 
Example  Problem  1 


DATA  FlLb  DtSCKlHllON  (47  ChAliACTEBS  MAX): 

INCLUDii  'CAT§PASS/FlLENAME$PAi>b '  ONLY,  NU  USLR  ID 
IF  YOU  VvlSri  lO  INPin  DATA  FROM  ItKMlNAi,  TYPE 
CAKRlAGfc,  KETURN. 

V 

ENTER  NEW  NAFiE  FOR  DATA  FII£  TO  BE  CREATED,  OR 
TYPE  CARRIAGE  RETURN  IF  YOU  DO  NOT  WISH  TO  SAVE  DATA. 
?  liUNJ 

ENTER  OUTPUT  FILE  NAME— CARRIAGE  RETURN  PRIMPS 
OUiPUl  AT  TERMINAL. 

?  OlUNj 

DATA  IS  ENIERED  IN  GROUPS  THAT  WILL  BE  DESCRIBED 
IN  THE  ORDER  OF  THEIR  ENTRY— ALONG  WITH  AN 
IDEWriFYlNG  SWBOL  Fa<  EACH  GROUP.  SEPARATE 
DATA  ON  EACH  LINE  BY  COMMAS  OR  AT  LEAST  ONE 
SPACE.  YOU  [AAY  REQUEST  SHORT  CUES  FCX<  IDEMITFYING 
SY^lbOL&  ONLY.  SYMBODd  ARE: 

*  =  TITLE  LINE 

#  =  NOLAL  POINT  DATA 
%  =  ELEMEMi  LATA 

+  =  NOLAL  POINT  LOAD 
$  =  ELEl'iEMi  U>'\D 
&  =  SPECIFICATION  LINE 


SEE  BRIEF  IXIR  DETAILED  INFORMATION  ON  INPUT  DATA 


SHORT  Oh  LONG  CUESi  (SH  OR  DO) 
?  LO 


GROUP  *  :  ENiER  TTTIE  FOR  THIS  PROBLEM  ON 
NEXI^  THREE  LINES.  (66  CHAR.  MAX.  PER  LINE) 

ENiER  SEND  IN  FIRST  LINE  IF  YOU  DO  NOT  WISH  TO 
CONTINUE. 

*  ?  NEWTUN  FOR  THE  DESIGN  OF  TUNNEIE 

*  ?  13  INCH  CONCRETE  LINER,  22  FT.  DIAMETER  TUNNEL 

*  ?  PARK  RIVER  DESIGN  (Na-llNAL  SUPPORT) 


y 

1 


(Continued) 
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Table  B1  (Continued) 


GROUP  &  : FIRST  SPECIFICATION  LINE 

ENTER  NO.  Of'  :OEAL  PDIMIS,  NO.  OF  EUiMENTS, 

ANL  NO.  OF  LOADING  CONUilTONS. 

&  ?  19  IB  1 


GROUP  #  ;  ENl'ER  NOEAL  POINT  LATA  AS  FOLDCMS: 
NODAL  Pi.  NO.,  REisTBAlNT  CODE  (3  DIGIT),  X-COOR, 
Y-COQR,  DISPL  IN  X-DIR,  DISPL  IN  Y-DIR 
ROTATION  ABCX/i  NODE  (RADIANS) . 

#  ?IlDIUt)k)k)i£) 

#  ?  2  id  2.m  id. 18  k)  Id  id 

«  V  3  0  3.95  0.70  000 

»  V  4  0  5.77  1.55  0  0  0 

#  ?  5  0  7.42  2.70  000 

#  ?  6  0  B.B4  4.12  000 

#  ?  7  0  9.99  5.77  000 

#  ?  8  0  10.85  7.59  0  0  0 

#  ?  9  0  11.37  9.54  000 

#  ?  10  0  11.54  il.54  000 

#  V  11  0  11.37  13.55  000 

#  ?  12  0  10.85  15.49  000 

#  ?  13  0  9.99  17.31  0  w  0 

«  ?  14  0  8.84  18.96  000 

#  ?  15  0  7.41  2w. 38  0  0  0 

#  ?  16  0  5.77  21.54  000 

#  ?  17  0  3.95  22.39  000 

#  ?  18  0  2.00  22.91  000 

#  ?  19  I0l  0.0  23.08  000 


GROUP  %  :  ENTER  ELEMENT  LATA  AS  FOLIOvS: 
EIEMENT  NO. ,NODfc:  NO.  1  END, NODE  NO.  J  END 
CROSS-SECT.  AREA, YOUNG'S  MODULUS 
FDN.  ROCK  MODULUS, MOMENl’  OF  INERTIA 
FI£XURAL  STIFF.  1  END, FLEXURAL  STIFF.  J  END, 
FI£XURAL  CARRY-OVER  FROM  I  END  TO  J  END 
%  ?  1  1  2  1.08  519119  1000.0  0.106  000 

%  ?  18  18  19  1.08  519119  1000.0  0.106  000 


ENTER  TITLE  FOR  LOADING  CONDITION  1 
*  ?  ROCK  LOAD  ONLY 


(Continued) 
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Table  B1  (Concluded) 


GROUP  &  : SPECIFICATION  LINE  FOR  LOADING  CONDITION  1 

B^rrE;R  no.  of  nodal  pi*,  load  line  eni'kies  to  polion, 

AND  A  ZERO  IF  YOU  WANT  ITERATIONS  PERFORMED 
TO  REMOVE  TENSION  IN  SPRINGS, OIHEFMISE, 

ENTER  A  -1  TO  EUMINATT;  ITERITIVE  PROCESS. 

«.  ?  li)  0 


GROUP  +  lENTER  NODAL  PT.  LCAD6  AS  FOLLOWS: 
NODAL  POINI  NO.  ,PC»CE  IN  X-DIR,PC»CE  IN  Y-DIK, 
MOMENT  ABOUT  NODE. 

+  ?  10  0  -0.09  0 

+  ?  11  0  -0.35  0 

+  ?  12  w  -0.69  0 

+  ?  13  0  -1.00  0 

+  ?  14  0  -1.29  0 

+  ?  15  0  -1.54  0 

+  ?  16  0  -1.74  0 

+  ?  17  0  -1.88  0 

+  ?  18  0  -1.97  0 

+  ?  19  0  -1.00  0 


GROUP  S  : ENTER  ELEMENT  LOAD  DATA  AS  FOLLOWS: 
ELEMENT  NO. , LOAD  TYPE  CODE, MAGNITUDE  OF  LOAD, 
LOAD  LOCATION, FORCE  DIRECTION  COLE 
a'YPE  ALL  zero's  TO  TERMINATE  LOAD  DATA) 

S  7  0  0  0  0  0 


DO  YOU  WANT  SHORT  OR  LONG  CUES  FOR  THIS  SECTION? 
?  SH 


?  SEND 


WOULD  YOU  LIKE  ANOTHER  RUN?  (YES  OR  NO) 
?  NO 


* 
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Table  B2 

Listing  of  Output  File  Created  from  Execution  Sequence 
For  Example  Problem  1  Shown  in  Table  B1 


LIST  ai‘UN3 


**  NEWl-UN  **  ANALYSIS  OF  CAST  IN  PLACE  TUNNELS 

^CMTUN  FOR  TilE  DESIO^  OF  TUNNEIS 

13  INCH  CONCRETE  LINER,  22  FT.  DIAMETER  TUNNEL 

PARK  RIVER  DESIGN  (NOMINAL  SUPPORT’) 

NUMBER  OF  NUDES  =  19 

NUMBER  OF  EUMENTS  =  18 
NUMBER  OF  LOADINGS  »  1 


(RESTRAINT'  CODE 

RESTRAINT' 

«1  SPECIFY  RESTRAINT, 

=8  SPECIFY  LOAD) 

SPECIFIED  DISPLACBMEMIS 

NODE 

CODE 

X-COOR 

Y-COOR 

X 

Y 

ROT 

1 

181 

0. 

8. 

0. 

8. 

0. 

2 

8 

2.00 

8.18 

8. 

0. 

0. 

3 

8 

3.95 

8.70 

0. 

0. 

0. 

4 

0 

5.77 

1.55 

0. 

8. 

8. 

6 

8 

7.42 

2.78 

0. 

0. 

0. 

6 

0 

8.84 

4.12 

0. 

8. 

0. 

7 

8 

9.99 

5.77 

8. 

0. 

0. 

8 

0 

10.85 

7.59 

0. 

8. 

0. 

9 

8 

11.37 

9.54 

0. 

0. 

8. 

lii 

8 

11.54 

11.54 

0. 

0. 

8. 

11 

8 

11.37 

13.55 

8. 

8. 

8. 

12 

0 

10.85 

15.49 

0. 

8. 

0. 

13 

8 

9.99 

17.31 

0. 

8. 

0. 

14 

0 

8.84 

18.96 

8. 

0. 

0. 

15 

8 

7.41 

28.38 

0. 

8. 

8. 

16 

8 

5.77 

21.54 

0. 

0. 

0. 

17 

8 

3.95 

22.39 

8. 

0. 

0. 

18 

0 

2.88 

22.91 

0. 

8. 

0. 

19 

181 

8. 

23.88 

0. 

0. 

8. 

(Continued)  (Sheet  1 
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Table  B2  (Continued) 


ELEMENT 

1 

J 

AREA 

EIASTIC-MOD 

SPRING  CONST 

1 

1 

2 

0.110000E 

01 

0.519119E 

06 

0.100000E 

04 

2 

2 

3 

0.110000E 

01 

0.519119E 

06 

0.100000E 

04 

3 

3 

4 

0. 110000E 

01 

0.519119E 

06 

0.100000E 

04 

4 

4 

5 

0.110000E 

01 

0.519119E 

06 

0. 100000E 

04 

5 

5 

6 

0.110000E 

01 

0.519119E 

06 

0.100000E 

04 

6 

6 

7 

0.110000E 

01 

0.519119E 

06 

0.100000E 

04 

7 

7 

8 

0.110000E 

01 

0.519119E 

06 

0.100000E 

04 

8 

8 

9 

0.110000E 

01 

0.519119E 

06 

0.100000E 

04 

9 

9 

10 

0. 110000E 

01 

0.519119E 

06 

0.100000E 

04 

Ik) 

10 

11 

0.110000E 

01 

0.519119E 

06 

0.100000E 

04 

11 

11 

12 

0.110000E 

01 

0.519119E 

06 

0.100000E 

04 

12 

12 

13 

0.110000E 

01 

0.519119E 

06 

0.100000E 

04 

13 

13 

14 

0.110000E 

01 

0.519119E 

06 

0.100000E 

04 

14 

14 

15 

0.110000E 

01 

0.519119E 

06 

0.100000E 

04 

15 

15 

16 

0.110000B 

01 

0.519119E 

06 

0. 100000E 

04 

16 

16 

17 

0.110000E 

01 

0.519119E 

06 

0.100000E 

04 

17 

17 

18 

0.110000E 

01 

0. 519119B 

06 

0.100000E 

04 

18 

18 

19 

0.110000B 

01 

0.519119E 

06 

0.100000C 

04 

E1£M£M' 

MOM  INERL'IA 

STIFF  IJ 

STIFF  JI 

COF  IJ 

COF  JI 

1 

0.100000E 

00 

4.000000 

4.000000 

0.500000 

0.500000 

2 

0.100000E 

00 

4.000000 

4.000000 

0.500000 

0.500000 

3 

0.100000E 

00 

4.000000 

4.000000 

0.500000 

0.500000 

4 

0.100000E 

00 

4.000000 

4.000000 

0.500000 

0.500000 

5 

0.100000E 

00 

4.000000 

4.000000 

0.500000 

0.500000 

6 

0.100000E 

00 

4.000000 

4.000000 

0.500000 

0.500000 

7 

0.100000E 

00 

4.000000 

4.000000 

0.500000 

0.500000 

8 

0.100000E 

00 

4.000000 

4.000000 

0.500000 

0.500000 

9 

0.100000E 

00 

4.000000 

4.000000 

0.500000 

0.500000 

10 

0. 100000E 

00 

4.000000 

4.000000 

0.500000 

0.500000 

11 

0.100000E 

00 

4.000000 

4.000000 

0.500000 

0.500000 

12 

0.100000E 

00 

4.000000 

4.000000 

0.500000 

0. 500000 

13 

0.100000E 

00 

4.000000 

4.000000 

0.500000 

0.500000 

14 

0.100000E 

00 

4.000000 

4.000000 

0.500000 

0.500000 

15 

0.100000E 

00 

4.000000 

4.000000 

0.500000 

0.500000 

16 

0.100000E 

00 

4.000000 

4.000000 

0.500000 

0.500000 

17 

0.100000E 

00 

4.000000 

4.000000 

0.500000 

0. 500000 

18 

0.100000E 

00 

4.000000 

4.000000 

0.500000 

0.500000 

(Continued) 
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Table  B2  (Continued) 


ROCK  LC^  ONLY 


NUMBER  OF  NOCAL  POINT  LOAD  CARC6  = 

SPECIFIED  NODAL  POINl'  LOADS 

NODE  FX  FY 

110  Id.  -Id.  09k) 

11  0.  -0.350 

12  0.  -0.690 

13  0.  -1.000 

14  0.  -1.290 

15  0.  -1.540 

16  0.  -1.740 

17  0.  -1.880 

18  0.  -1.970 

19  0.  -1.000 


ELEMENT  LOADS 

LOAD  TYPE  *  *1  CONCENT’RATED  LOAD 

=2  UNIFORM  DISTRIBUTED  LOAD 
“3  CCNCENl'RATED  MOMENT 
=4  TRIANGULAR  DISTRIBUTED  LOAD 

LOCATION  CC®E-  *LQAD  TYPES  1  AND  3: 

DISTANCE  FROM  I  NODE 
*LQAD  TYPE  2:  =0 
*LQAD  TYPE  4; 

=0  NODE  I  WEIGHTED,  >1  NODE  J  WEIGHTED 


DIRECTION  CODE- 


*LQAD  TYPES  1-4: 
X»2,  Y=1 


element  node  LEAD  MASilTUDE  LOCATION  DIRECTION 
I  J  TYPE  OF  LOAD  CODE  CODE 


(Continued) 
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Table  B2  (Continued) 


Eimh^lS  ON  AN  £1ASTIC  FOUNDATION 


£I£MENr  SPRING  PRESSURE 

1  J 


1 

8.111415E 

01 

0.113347E 

01 

2 

8.189222E 

01 

0.115354E 

01 

3 

8.186871E 

01 

0.117264E 

01 

4 

0.104883E 

01 

0. 119862E 

01 

5 

8.102877E 

01 

0.122638E 

01 

6 

0.101346E 

01 

0.126574E 

01 

7 

0.102393E 

01 

0. 133689B 

01 

8 

0.103027E 

01 

0.139396E 

01 

9 

0.ie4996E 

01 

0.143565E 

01 

18 

0.10 607 5E 

01 

0.141759E 

01 

11 

0.978173E 

00 

0.120488E 

01 

12 

0.717539E 

00 

0.691237E 

00 

13 

0.209308E 

00 

0. 

14 

0. 

0. 

15 

0. 

0. 

16 

0. 

0. 

17 

0. 

0. 

18 

0. 

0. 

NUMBER  OP 

ITERATIONS  * 

3 
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Table  B2  (Continued) 


El£HENr 

ZONE 

LIMITS 

OP  INTEGRATION  (HAX=20) 

1 

ODD 

0. 

9 

2.008, 

2 

ODD 

0. 

9 

2.018, 

3 

ODD 

0. 

9 

2.009, 

4 

ODD 

0. 

9 

2.011, 

5 

ODD 

0. 

9 

2.008, 

6 

ODD 

0. 

9 

2.011, 

7 

ODD 

0. 

9 

2.013, 

8 

ODD 

0. 

9 

2.018, 

9 

ODD 

0. 

9 

2.007, 

10 

ODD 

0. 

9 

2.017, 

11 

ODD 

0. 

9 

2.008, 

12 

ODD 

0. 

9 

2.013, 

13 

ODD 

0. 

9 

1.320,  2.011, 

14 

EVEN 

0. 

9 

2.015, 

15  EVEN 

0. 

9 

2.009, 

16  EVEN 

0. 

9 

2.009, 

17  EVEN 

0. 

9 

2.018, 

18  EVEN 

0. 

9 

2.007, 

NOCAL  POINT  DISPLACEMENIS 


•DE 

X-DISPL 

Y-^JISPL 

ROT-DISPL 

1 

0. 

-0.11186E-02 

0. 

2 

-0.'<t3378E-04 

-0.11420E-02 

-0.19741E-04 

3 

-0.71554E-04 

HS.12129E-02 

-0.40994E-04 

4 

-0.68622E-04 

-0.13263E-02 

-0.62894E-04 

5 

-0.20191EH04 

-0.14751E-02 

-0.86186E-04 

6 

0.87515E-04 

-0.16469E-02 

-0.11124E-03 

7 

0.26856E-03 

-0.18283E-02 

-0.14056E-03 

8 

0.53212E-03 

-0.20031E-02 

-0.16933E-03 

9 

0.87177E-03 

-0.214098-02 

-0.18931E-03 

10 

0.12522E-02 

-0.22190E-02 

-0.19088E-03 

11 

0.16116E-02 

-0.22345E-02 

-0.15472E-03 

12 

0.18422E-02 

-0.22192E-02 

-0.60764E-04 

13 

0.18370E-02 

-0.22696B-02 

0.95663E-04 

14 

0.15440E-02 

-0.25233E-02 

0.28811E-03 

15 

0. 10427E-02 

-0.30804E-02 

0.44065E-03 

16 

0.51758E-03 

-0.38791E-02 

0.49049E-03 

17 

0.14809E-03 

-0.47374E-02 

0.41897E-03 

18 

-0.26211E-05 

-0.54011E-02 

0.23998E-03 

19 

0. 

-0.56487E-02 

0. 
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Table  B2  (Concluded) 


END  ACTIONS  WITH  NOEAL  POINT  LCAD6 


El^MENl  AXIAL  I 

AXIAL  J 

1  12.881) 

-12.888 

2  12.878 
-12.878 

3  12.879 
-12.879 

4  12.879 
-12.879 

5  12.888 
-12.888 

6  12.876 
-12.876 

7  12.883 
-12.883 

S  12.912 

-12.912 

9  12.963 

-12.963 

18  12.964 

-12.964 

11  12.768 
-12.768 

12  12.287 
-12.287 

13  11.529 
-11.529 

14  18.446 
-18.446 

15  9.233 
-9.233 

16  8.875 
-8.075 

17  7.188 
-7.188 

18  6. 788 
-6.788 


SHEAR  I 

MOMENT  I 

SHEAR  J 

MOMENT'  J 

-1.159 

0.097 

-1.891 

-8.172 

-1.115 

0.172 

-1.144 

-0.163 

-1.132 

0.163 

-1.111 

-0.219 

-1.106 

0.219 

-1.146 

-0.229 

-1.136 

0.229 

-1.120 

-0.312 

-1.162 

0.312 

-1.121 

-0.438 

-1.839 

0.438 

-1.328 

-0.253 

-1.011 

0.253 

-1.428 

0.045 

-0.852 

-0.045 

-1.642 

8.708 

-0.556 

-0.708 

-1.956 

1.998 

-0.274 

-1.998 

-1.950 

3.605 

-0.866 

-3.605 

-1.484 

4.968 

-0.089 

-4.968 

-0.067 

5.823 

-1.085 

-5.823 

1.085 

2.836 

-1.541 

-2.836 

1.541 

-0.260 

-1.581 

0.260 

1.581 

-3.437 

-1.157 

3.437 

1.157 

-5.771 

-0.433 

5.771 

8.433 

-6.641 

***  $END  *** 
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Example  Problem  2 

5.  The  tunnel  section  used  in  this  example  is  the  same  as  that 
used  in  example  1.  The  four  loading  cases  that  are  used  are  delineated 
in  the  loading  diagram  shown  in  Figure  B4*  These  loading  cases  include 
the  weight  of  the  lining  material  which  was  not  considered  in  example  1. 

6.  Table  B3  shows  the  execution  of  the  program  from  a  data  file. 
The  data  file  needed  to  run  this  problem  is  shown  in  Table  B4.  The 
listing  of  the  output  file  is  given  in  Table  B5. 

7.  Figure  B5  gives  the  formula  used  in  calculating  the  nodal 
loads  due  to  the  rock  load.  Figure  B6  shows  sample  calculations  to 
obtain  the  nodal  load  for  node  16  for  load  case  1.  Figure  B7  gives  the 
formulas  for  calculating  nodal  loads  due  to  water  loading.  Figures  B8 
and  B9  show  sample  calculations  for  nodal  loads  at  node  16  for  load 
cases  2  and  3,  respectively.  Figure  BIO  shows  a  plot  of  the  moment  in 
the  concrete  lining  for  each  load  case  of  example  2. 
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A  =  108  ! 

4  1 

I  =  0.106  ft  I 

R  =  1,000  kef  ( 


LOAD 

CASE 

DESCRIPTION 

1 

ROCK  LOAD  +  WEIGHT  OF  LINING 

2 

ROCK  LOAD  +  WEIGHT  OF  LINING  +  30-ft  NET  INTERNAL  HEAD 

3 

ROCK  LOAD  +  WEIGHT  OF  LINING  +  85-ft  NET  EXTERNAL  HEAD 

4 

ROCK  LOAD  +  WEIGHT  OF  LINING  +  185-ft  NET  EXTERNAL  HEAD 

NOTE:  COMPUTER  MODEL  AND  ROCK  MODULUS  FOR  THIS  EXAMPLE 
ARE  THE  SAME  AS  THOSE  FOR  EXAMPLE  1  (SEE  FIGURES 
82  AND  83). 


Figure  B4.  Physical  problem  treated  in  example  2 
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185- ft  EXTERNAL 


Table  B3 

Execution  of  NEVnUN  from  the  Data  File 
For  Example  Problem  2  Shovm  in  Table  B4 


DATA  FILE.  DEiSCklPIlUN  (47  CriARACi’EftS  MAX) ; 

INCLUDE  'CA'r$FASS/FILENAME$PASS '  ONLY,  NO  USER  ID 
IF  YOU  WISH  TO  INPUI'  DATA  FKCM  TEFMINAL,  TYPE 
CARRIAGE  RETURN. 

?  TUNIN 

ENTER  OUl’PUr  F1L£  NAME— CARRIAGE  RETURN  PRINTS 
OUrPUI'  AT  TERMINAL. 

?  TUNOUr 


WOULD  YOU  LIKE  ANOl'HER  RUN?  (YES  OR  NO) 
?  NO 


lOuO  HEi.'TLih  Fi'i^MFLE  FPDFLE^l 
lOlU  13  IhCH  CONfFETE  LIhFP  cc:  FT  DIRf-T  Ti.ihHEL  (title  lines) 
Ovh  FIVEF  rEIIi.-fi  'riGEMFL  FCCK  EUPFnFT' 


_  _  (specification  line) 


load  case  sped 


14*:<0  1  0^-1. 

l^Mij  i  -1  . 

1510  3*0. - ; . 
I5v0 


^roi4CT^«,M.giff701 


1580  10« 1.881-. 
1580  111  1.871-. 


(nodal  load  data) 


(Continued) 
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load  case  specification  element  line 


load  case  title  line 


load  case  specification  node  line 


BnMWM^igBBiBgiTrfMigaroni 


load  case  title  line) 


load  case  specification  node  line 


A_..  _|L 


Table  B5 

Listing  of  Output  Results  for  Example  2 


1 


**  NEMTLN  **  ANALYSIS  OF  CAST  IN  PLACE  TIWNEIS 

NENTUN  EXAMPLE  PRCBIEM 

13  INCH  CONCRETE  LINER  22  FT  DIAM  TUNNEL 

PARK  RIVER  DESIGN  (NC«MAL  ROCK  SUPPORT) 

NUMBER  OF  NODES  =  19 

NUMBER  OF  EIEMENTS  «  18 
NUMBER  OF  LOADINGS  »  4 


(RESTRAINT’  CODE  =1  SPECIFY  RESTRAINT',  *0  SPECIFY  LOAD) 


RESTRAINT' 


DE 

CODE 

X-COC» 

Y-COQR 

1 

101 

0. 

0. 

2 

0 

2.00 

0.18 

3 

0 

3.95 

0.70 

4 

0 

5.77 

1.55 

5 

0 

7.42 

2.70 

6 

0 

8.84 

4.12 

7 

0 

10.00 

5.77 

8 

0 

10.85 

7.59 

9 

0 

11.37 

9.53 

10 

0 

11.54 

11.54 

11 

0 

11.37 

13.55 

12 

0 

10.85 

15.49 

13 

0 

10.00 

17.31 

14 

0 

8.84 

18.96 

15 

0 

7.42 

20.36 

16 

0 

5.77 

21.54 

17 

0 

3.95 

22.39 

18 

0 

2.00 

22.91 

19 

101 

0. 

23.08 

SPECIFIED  DISPLACBtENTS 


X  Y  ROT 

0.  0.  0. 

0.  0.  ). 

0.  0.  0. 

0.  0.  0. 

0.  0.  0. 

0.  0.  0. 

0.  0.  0. 

0.  0.  0. 

0.  0.  0. 

0.  0.  0. 

0.  0.  0. 

0.  0.  0. 

0a  0«  0* 

0.  0.  0. 

0.  0.  0. 

0.  0.  0. 

0.  0.  0. 

0.  0.  0. 

0.  0.  0. 
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Table  B5  (Continued) 


E1£M£M' 

I 

J 

A8LA 

EIASTIC-MOD 

SPRING  CONST 

1 

1 

2 

8.188888E 

01 

0.519119E 

06 

0.100000E 

04 

2 

2 

3 

i3.l88888E 

01 

0.519119E 

06 

0.100000B 

04 

3 

3 

4 

8.108808E 

01 

0.519119E 

06 

0.100000E 

04 

4 

4 

5 

id.immm 

01 

0.51911Sffi 

06 

0.100000E 

04 

5 

5 

6 

8.188008E 

01 

0.519119E 

06 

0.100000E 

04 

6 

6 

7 

Id.  108800E 

01 

0.519119E 

06 

0.100000E 

04 

7 

7 

8 

0.108800E 

01 

0.519119E 

06 

0.100000E 

04 

8 

8 

9 

0.108000E 

01 

0.519119E 

06 

0.100000E 

04 

9 

9 

18 

0.10b000E 

01 

0.519119E 

06 

0.100000E 

04 

18 

18 

11 

0.108000E 

01 

0.519119E 

06 

0. 100000E 

04 

11 

11 

12 

0.108000E 

01 

0.519119E 

06 

0.100000E 

04 

12 

12 

13 

0.108000E 

01 

0.519119E 

06 

0. 100000E 

04 

13 

13 

14 

0.108000E 

01 

0. 519119E 

06 

0.100000E 

04 

14 

14 

15 

0.108000E 

01 

0.5i9119E 

06 

0.100000E 

04 

15 

15 

16 

0.108000E 

01 

0.519119E 

06 

0. 100000E 

04 

16 

16 

17 

0.108000E 

01 

0.519119E 

06 

0.100000E 

04 

17 

17 

18 

0.l080b0E 

01 

e.519119E 

06 

0.100000E 

04 

18 

18 

19 

0.108000E 

01 

0.519119E 

06 

0.100000E 

04 

EUBMEMl' 

MCM  INEKIIA 

STIFF  IJ 

STIFF  JI 

CCF  IJ 

COF  JI 

1 

0.110000E 

00 

4.000000 

4.000000 

0.500000 

0.500000 

2 

0.110000E 

00 

4.000000 

4.000000 

0.500000 

8.508008 

3 

0.110000E 

00 

4.000000 

4.000000 

0.500000 

0.500000 

4 

0.110000E 

00 

4.000000 

4.000000 

0.500000 

0.500000 

5 

0.110000E 

00 

4.000000 

4.000000 

0.500000 

0.500000 

6 

0.110000E 

00 

4.000000 

4.000000 

0.500000 

0.580000 

7 

0.110000E 

00 

4.000000 

4.000000 

0.500000 

0.580000 

8 

0.110000E 

00 

4.000000 

4.000000 

0.500000 

0. 500000 

9 

0. 110000E 

00 

4.000000 

4.000000 

0.500000 

8.500000 

10 

0.110000E 

00 

4.000000 

4.000000 

0.500000 

0.500000 

11 

0. 110000E 

00 

4.000000 

4.000000 

8.500000 

0.500000 

12 

0.110000E 

00 

4.000000 

4.000000 

0.508000 

0.500000 

13 

0.110000E 

00 

4.000000 

4.000000 

0.500000 

0.508800 

14 

0.110000E 

00 

4.000000 

4.000000 

0.580000 

0.500000 

15 

0.110000E 

00 

4.000000 

4.000000 

0.500000 

8.500008 

16 

0.110000E 

00 

4.000000 

4.000000 

0. 500000 

0.500000 

17 

0.110000E 

00 

4.000000 

4.000000 

0.508000 

0.500000 

18 

0.110000E 

00 

4.000000 

4.000000 

0.500000 

0.500000 

LC-1  ROCK  LOAD  0.275r/PT 
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Table  B5  (Continued) 


NIMBEK  CF  NOCAL  POINT  1£AD  CARD6  »  19 


SPBCIFIED  NOIAL  POINT  1£AI]6 


NODE 

FX 

FY 

MOM 

1 

0. 

-k).16tj 

0. 

2 

0. 

-8.328 

0. 

3 

0. 

-8.320 

0. 

4 

0. 

-0.320 

0. 

5 

0. 

-0.320 

0. 

6 

0. 

-0.320 

0. 

7 

0. 

-0.320 

0. 

8 

0. 

-0.320 

0. 

9 

0. 

-0.320 

0. 

It) 

0. 

-0.360 

0. 

11 

0. 

-0.630 

0. 

12 

0. 

-0.710 

0. 

13 

0» 

-0.900 

0. 

14 

0. 

-1.070 

0. 

15 

0. 

-1.210 

0. 

lb 

0. 

-1.320 

0. 

17 

0. 

-1.410 

0. 

18 

0. 

-1.470 

0. 

19 

0. 

-0.820 

0. 

ELEMEl^r^  LOADS 

LOAD  TYPE  *  =1 

CONCENTRATED  LOAD 

=2 

UNIFa^% 

DISTRIBUTED  LOAD 

=3  CONCENTRATED  MOMENT 
=4  TRIANGUIAR  DISTRIBUTED  LOAD 

LOCATION  CODE-  *LOAD  TYPES  1  AND  3: 

DISTANCE  PRCM  1  NODE 
♦LOAD  TYPE  2;  =0 
♦LOAD  TYPE  4: 

=0  NODE  1  MEIGHTEOr  NODE  J  WEIOITED 

DIRECTION  CODE-  *LOAD  TYPES  1-4; 

X=2,  Y=1 
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Table  B5  (Continued) 


EI£MEMl‘  NODE  LOAD  MAGNITUDE  UXATIOl  DIRECTION 
I  J  TYPE  OF  LOAD  CODE  CODE 


EIEMEMTS  ON  AN  EIASTIC  FOUNDATION 

ElEHENr  SPRING  PRESSURE 

I  J 


1 

0.121958E 

01 

0.123552E 

01 

2 

0.119140E 

01 

0. 124252E 

01 

3 

8.115213E 

01 

0.123929E 

01 

4 

Id.  11881 7E 

01 

0.123479E 

01 

5 

8.185648E 

01 

0. 12251 0E 

01 

6 

8.180927E 

01 

0.122624E 

01 

7 

0.965842E 

00 

0.122936E 

01 

8 

0.928718E 

00 

0.123352E 

01 

9 

8.887968E 

00 

0.121106E 

01 

lid 

8.845478E 

00 

0.113649E 

01 

11 

8.715821E 

00 

0.888689E 

00 

12 

0.442099E 

00 

0.398478E 

00 

13 

0. 

0. 

14 

0. 

0. 

15 

0. 

0. 

16 

0. 

0. 

17 

0. 

0. 

18 

0. 

0. 

NUMBER  OF 

ITERATIONS 

^  * 

3 

(Continued) 
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Table  B5  (Continued) 


E1£MCNT  ZONE  LIMIIS 

1  ODD  , 

2  01X3  0.  , 

3  ODD  Id.  , 

4  ODD  a.  , 

5  ODD  0.  , 

6  ODD  0.  , 

7  ODD  0.  , 

8  01X3  0.  , 

9  ODD  0.  , 

10  ODD  0.  , 

11  ODD  0.  , 

12  ODD  0.  , 

13  EVEN  0.  , 

14  EVEN  0.  , 

15  EVEN  0.  , 

16  EVEN  0.  , 

17  EVEN  0.  , 

18  EVEN  0.  , 


NODAL  fOINT  DISPLACEMENTS 
NODE  X-DISPL 


1 

0. 

-0 

2 

-0.42323E-04 

-0 

3 

-0.71616E-04 

-0 

4 

-0.74014E-04 

-0 

5 

-0.36755E-04 

-0 

6 

0.52815E-04 

-0 

7 

0.20653E-03 

-0 

8 

0.42900E-03 

-0 

9 

0.71326E-03 

-0 

10 

0.10319E-02 

-0 

11 

0.13251E-02 

-0 

12 

0.15021E-02 

-0 

13 

0.14770E~02 

-0 

14 

0.12196E-02 

-0 

15 

0.81817E-03 

-0 

16 

0.39953E-03 

-0 

17 

0.11437E-03 

-0 

18 

-0.16338E-05 

-0 

19 

0. 

-0 

OP  INTEGRATION  (MAX-20) 

2.008, 

2.018, 

2.009, 

2.011, 

2.008, 

2.017, 

2.009, 

2.008, 

2.017, 

2.017, 

2.008, 

2.009, 

2.017, 

1.994, 

2.029, 

2.009, 

2.018, 

2.007, 


Y-DISPL  ROT-DISPL 

12244E-02  0. 

12443E-02  -0 . 16407E-04 

13050E-02  -0.34280E-04 

14024E-02  -0.52947E-04 

15307E-02  -0.73219E-04 

16797E-02  -0.95472E-04 

18383E-02  -0.11990E-03 

19866E-02  -0.14292E-03 

21034E-02  -0.15886E-03 

21690E-02  -0.15786E-03 

21818E-02  -0.12368E-03 

21716E-02  -0.40906E-04 

22209E-02  0.90766E-04 

24406E-02  0.24121E-03 

28876E-02  0.34927E-03 

35160E-02  0.38007E-03 

41775E-02  0.32187E-03 

46871E-02  0.18443E-03 

48777E-02  0. 
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Table  B5  (Continued) 


1 


END  ACTIONS  WITH  NODAL  POINT  LOADS 


£1£MEN1' 

1 

2 

3 

4 

5 

6 

7 

8 
9 

lid 

11 

12 

13 

14 

15 

16 

17 

18 


AXIAL  I 
AXIAL  J 

12.267 
-12.267 
12. 2i)9 
-12.269 
12. 106 
-12. 100 
11.941 
-11.941 
11.736 
-11.736 
11.491 
-11.491 
11.225 
-11.225 
10.950 
-10.950 
10.687 
-10.687 
10.417 
-10.417 
10.033 
-10.033 
9.509 
-9. 509 
8.797 
-8.797 
7.872 
-7.872 
6.916 
-6.916 
6.017 
-6.017 
5.337 
-5.337 
4.965 
-4.965 


SHEAR  I 
SHEAR  J 

-1.265 

-1.195 

-1.217 

-1.233 

-1.218 

-1.178 

-1.171 

-1.179 

-1.165 

-1.118 

-1.098 

-1.148 

-1.015 

-1.181 

-0.875 

-1.291 

-0.690 

-1.428 

-0.357 

-1.653 

-0.075 

-1.564 

0.093 

-0.982 

-0.179 

0.179 

-0.987 

0.987 

-1.225 

1.225 

-1.273 

1.273 

-0.947 

0.947 

-0.401 

0.401 


MOMENT  1 
MOMENT  J 

0.017 

-0.093 

0.093 

-0.095 

0.095 

-0.165 


0.199 

-0.303 

0.303 

-0.326 

0.326 

-0.249 

0.249 

0.067 

-0.067 

0.702 

-0.702 

1.910 

-1.910 

3.347 

-3.347 

4.440 

-4.440 

4.079 

-4.079 

2.110 

-2.110 

-0.376 

0.376 

-2.934 

2.934 

-4.844 

4.844 

-5.649 


LC:-2  ROCK  LOAD  PLUS  INTERNAL  HEAD 
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Table  B5  (Continued) 


NljMBER  OP  NOIAL  POIMT  LCAD  CARDS 
SPECIFIED  NODAL  POINT  LOADS 


>OE 

FX 

FY 

1 

0. 

-1.000 

2 

0.290 

-1.970 

3 

0.580 

-1.910 

4 

0.860 

-1.810 

5 

1.120 

-1.660 

6 

1.360 

-1.460 

7 

1.580 

-1.230 

8 

1.750 

-0.960 

9 

1.880 

-0.650 

10 

1.960 

-0.360 

11 

1.970 

-0.180 

12 

1.930 

-0.010 

13 

1.810 

0.150 

14 

1.630 

0.300 

15 

1.390 

0.450 

16 

1.100 

0.580 

17 

0.790 

0.670 

18 

0.390 

0.730 

19 

0. 

0.290 

EL£MEl«r  LJADS 

LOAD  TYPE  *  =1  CaCEMfHATED  LOAD 

=2  UNIFORM  DISTRIBWED  LCAD 
*3  CONCEMl'RATED  MOMENT 
M  TRIANGUIAR  DISTRIBUTED  LOAD 

LOCATION  CODE-  *LOAD  TYPES  1  AND  3; 

DISTANCE  FROM  I  NODE 
♦LOAD  TYPE  2:  -0 
♦LOAD  TYPE  4: 

^  NODE  I  WEIQITED,  *1  NODE  J  WEICSITED 

DIRECTION  CODE-  *LOAD  TYPES  1-4: 

X-2,  Y=1 


ELEMENT'  NODE  LOAD  MAGNITUDE  LXATION  DIRECTION 
I  J  TYPE  OF  LOAD  CODE  CODE 
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Table  B5  (Continued) 


ELEMENTS  CM  AN  ELASTIC  PCXJNDATICN 


ELEMENT  SPRING  PRESSURE 

I  J 


1 

0.993190E 

00 

0.100527E 

01 

2 

0.975429E 

00 

0.101257E 

01 

3 

0.951448E 

00 

0.101391E 

01 

4 

0.925156E 

00 

0.101422E 

01 

5 

0.893329E 

00 

0.101125E 

01 

6 

0.864796E 

00 

0.101551E 

01 

7 

0.838110E 

00 

0.102141E 

01 

8 

0.816958E 

00 

0.102888E 

01 

9 

0.793574E 

00 

0. 101993E 

01 

10 

0.770731E 

00 

0.978370E 

00 

11 

0.690692E 

00 

0.827191E 

00 

12 

0. 52157 3E 

00 

0.519707E 

00 

13 

0.198525E 

00 

0. 

14 

0. 

0. 

15 

0. 

0. 

16 

0. 

0. 

17 

0. 

0. 

18 

0. 

0. 

NUMBER  OF  ITERATIONS  »  3 
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Table  B5  (Continued) 


i 


EUIMEM' 

ZONE 

LLMllS  OF 

INTEGRATION  (MAX- 

1 

ODD 

0. 

r 

2.008, 

2 

ODD 

0. 

9 

2.018, 

3 

ODD 

0. 

9 

2.009, 

4 

ODD 

0. 

9 

2.011, 

5 

OW) 

0. 

9 

2.008, 

6 

ODD 

0. 

9 

2.017, 

7 

ODD 

0. 

9 

2.009, 

8 

ODD 

0. 

9 

2.008, 

9 

ODD 

0. 

9 

2.017, 

10 

ODD 

0. 

9 

2.017, 

11 

ODD 

0. 

9 

2.008, 

12 

ODD 

0. 

9 

2.009, 

13 

ODD 

0. 

9 

1.979,  2.017, 

14 

EVEN 

0. 

9 

1.994, 

15  EVEN 

0. 

9 

2.029, 

16 

EVEN 

0. 

9 

2.009, 

17  EVEN 

0. 

9 

2.018, 

18 

EVEN 

0. 

9 

2.007, 

f 


NOCAL 

POINT  DISPLACEMENTS 

NODE 

X-OISPL 

V-DISPL 

RCa’-DlSH. 

1 

0. 

-0.99720E-03 

0. 

2 

0.10119E-05 

-0.10892E-02 

-0.12118E-04 

3 

0.10714E-04 

-0.10451E-02 

-0.24759E-04 

4 

0.37617E-04 

-0.11015E-02 

-0.37559E-04 

5 

0.89437E-04 

-0.11739E-02 

H8.51300E-04 

6 

0.17411E-03 

-0.12560E-02 

-0.66507E-04 

7 

0.29896E-03 

-0.13405E-02 

-0.83154E-04 

8 

0.46659E-03 

-0. 14147E-02 

H8.99153E-04 

9 

0.67250E-03 

-0.14651E-02 

-0.11080E-03 

10 

0.89853E-03 

-0.14785E-02 

-0.1111ffi:-03 

11 

0.11049E-02 

-0.14542E-02 

—0 . 90688E“04 

12 

0.12345E-02 

-0.14108E-02 

-0.39669E-04 

13 

0.12285E-02 

-0.14022E-02 

0.46817E-04 

14 

0. 10538E-02 

-0.15092E-02 

0. 15670E-03 

15 

0.75397E-03 

-0.17907E-02 

0.24879BH93 

16 

0.41620E-03 

-0.22297E-02 

0.28440E-03 

17 

0.16552E-03 

-0.27138E-02 

0.24574E-03 

18 

0.38269E-04 

-0.30947E-02 

0.1421«:-03 

19 

0. 

-0.32385E-02 

0. 
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Table  B5  (Continued) 


# 


END  ACniONS 

WITH  NODAL  POINT 

ir:Ai£ 

ElEMENT 

AXIAL  1 

SHEAR  I 

MGMENI  1 

AXIAL  J 

SHEAR  J 

MOMENT  J 

i 

0.020 

-1.006 

0.003 

-0.020 

-0.997 

-0.016 

2 

-0.038 

-0.999 

0.016 

0.038 

-1.002 

-0.026 

3 

-0.145 

-0.983 

0.026 

0.145 

-0.987 

-0.043 

4 

-0.303 

-0.979 

0.043 

0.303 

-0.966 

-0.086 

5 

-0.510 

-0.961 

0.086 

0.510 

-0.946 

-0.141 

6 

-0.753 

-0.933 

0.141 

0.753 

-0.958 

-0.167 

7 

-1.019 

-0.877 

0.167 

1.019 

-0.985 

-0.120 

B 

-1.306 

-0.791 

0.120 

1.306 

-1.058 

0.077 

9 

-1.588 

-0.656 

-0.077 

1.588 

-1.173 

0.522 

10 

-1.892 

-0.499 

-0.522 

1.892 

-1.272 

1.230 

11 

-2.322 

-0.270 

-1.230 

2.322 

-1.271 

2.188 

12 

-2.891 

-0.088 

-2.188 

2.891 

-0.986 

3.091 

13 

-3.592 

-0.090 

-3.091 

3.592 

-0.142 

3.212 

14 

-4.459 

-0.577 

-3.212 

4.459 

0.577 

2.062 

15 

-5.363 

-1.045 

-2.062 

5.363 

1.045 

-0.057 

16 

-6.215 

-1.037 

0.057 

6.215 

1.037 

-2.141 

17 

-6.891 

-0.782 

2.141 

6.891 

0.782 

-3.719 

18 

-7.249 

-0.325 

3.719 

7.249 

0.325 

-4.372 

LC-3  RXK  LOAD  PLUS  EXTERNAL  HEAD 
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Table  B5  (Continued) 


NUMBER  CJF 

NODAL  POINT 

II 

O 

Q 

19 

PECiFIED 

NODAL  POINT 

LOADS 

NODE 

FX 

FY 

MOM 

i 

0. 

6.090 

0. 

2 

-2.170 

11.990 

0. 

3 

-4.270 

11.420 

0. 

4 

-6.240 

10.490 

0. 

-8.00kJ 

9.220 

0. 

0 

-9.510 

7.660 

0. 

1 

-10.710 

5.860 

0. 

b 

-11.590 

3.900 

0. 

9 

-12.090 

1.810 

0. 

10 

-12.240 

-0.360 

0. 

11 

-12.000 

-2.650 

0. 

12 

-11.410 

-4.860 

0. 

13 

-10. 480 

-6.950 

0. 

14 

-9.240 

-8.820 

0. 

15 

-7.630 

-10.430 

0. 

16 

-6.000 

-11.710 

0. 

17 

-4.100 

-12.670 

0. 

18 

-2.080 

-13.250 

0. 

19 

0. 

-6.850 

0. 

EI£MEN1'  IDACb 

liJAD  TYPE  *  =1  CONCENTRATED  LOAD 

=2  UNIFORM  DISTRIBUTED  LOAD 
=3  CONCENl'RATED  MOMENT 
=4  TRIANGULAR  DISTRIBUTED  LOAD 

LOCATION  CODE-  *LQAD  TYPES  1  AND  3: 

DISTANCE  FRCM  1  NODE 
*LCAD  TYPE  2:  =0 
♦LOAD  TYPE  4: 

=0  NODE  1  WEIGHTED,  =1  NODE  J  WEIGHTED 

DIRECTION  code-  ♦LOAD  TYPES  1-4; 

X=2,  Y=i 

ELEMENT  NODE  LOAD  MAGNITUDE  LXATION  DIRECTION 
1  J  TYPE  OF  LOAD  CODE  CODE 
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Table  B5  (Continued) 


El£MENlb  CM  AN  ELASTIC  FOUNDATION 


EIEMENI  SPRING  PRESSURE 

I  J 


1 

k).899391E 

00 

0.928743E 

00 

2 

8.852459E 

00 

0.955864E 

00 

3 

^.798761E 

00 

0. 97992 5E 

00 

4 

0.751198E 

00 

0.102243E 

01 

5 

l().709362E 

00 

0.107185E 

01 

6 

0.690016E 

00 

0.115451E 

01 

7 

0.687319E 

00 

0.123777E 

01 

8 

0.693635E 

00 

0.129725E 

01 

9 

k).664313E 

00 

0.125327E 

01 

lu 

13.575781E 

00 

0.103303E 

01 

11 

0.24 581 7E 

00 

0.385992E 

00 

12 

0. 

13 

0. 

0. 

14 

0. 

0. 

15 

0. 

0. 

16 

0. 

0. 

17 

0. 

0. 

18 

0. 

0. 

NUMBER  CF  ITERATIC»iS  =  4 
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Table  B5  (Continued) 


EUiMENl 

ZONE 

LIMITS 

OF  IWTEC 

1 

ODD 

0. 

f 

2.008 

2 

ODD 

0. 

9 

2.018, 

3 

ODD 

0. 

9 

2.009 

4 

ODD 

0. 

9 

2.011 

5 

ODD 

id. 

9 

2.008 

6 

ODD 

id. 

9 

2.017 

7 

ODD 

d. 

9 

2.009 

8 

ODD 

d. 

9 

2.008 

9 

ODD 

0. 

9 

2.017 

10 

ODD 

0. 

9 

2.017 

11 

ODD 

d. 

9 

2.008^ 

12 

EVEN 

id. 

9 

2. 009 

13 

EVEN 

0. 

9 

2.017 

14 

EVEN 

0. 

9 

1.994 

15  EVEN 

d. 

9 

2.029 

16  EVEN 

id. 

9 

2.009 

17 

EVEN 

d. 

9 

2.018 

Ifa 

EVEN 

d. 

9 

2.007 

NODAL  roiNTi  DISPLACEMENTS 


I 


DE 

X-DISPL 

Y-DISPL 

ROT-DISH. 

1 

0. 

-0.90303E-03 

0. 

2 

-0.28443E-03 

-0.95810E-03 

-0.31835E-04 

3 

-0.53745E-03 

-0.11326E-02 

-0.70219E-04 

4 

-0.72144E-03 

-0.14185E-02 

-0.11197E-03 

5 

-0.80215E-03 

-0.18053E-02 

-0.15732E-03 

6 

-0.74823E-03 

-0.22641E-02 

-0.20527E-03 

7 

-0. 53306E-03 

-0.27656E-02 

-0.25415E-03 

8 

-0.15471E-03 

-0.32563E-02 

-0.29161E-03 

9 

0.35488E-03 

-0.36866E-02 

-0.30357E-03 

10 

0.91779E-03 

-0.40195E-02 

-0.27247E-03 

11 

0. 13973E-02 

-0.42636E-02 

-0.16394E-03 

12 

0. 16056E-02 

-0.44992E-02 

0.35989E-04 

13 

0.14395E-02 

-0.48854E-02 

0.27726E-03 

14 

0.96017E'03 

-0.55618E-02 

0.48513E-03 

15 

0.38830E-03 

-0.65275E-02 

0. 58370E-03 

16 

-0.81659E-04 

-0.76513E-02 

0.57469E-03 

17 

-0.29447E-03 

-0.87340E-02 

0.47300E-03 

18 

-0.23661E-03 

-0.9540SE-02 

0.27005E-03 

19 

0. 

-0.98411E-02 

0. 

(Continued) 
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Table  B5  (Continued) 


END  ACTIONS  WITH  NOCAL  POINT  IJDAE6 


ELEMENT 

AXIAL  I 

SHEAR  I 

MCMENI'  1 

AXIAL  J 

SHEAR  J 

MOMENT  J 

1 

80.469 

-1.128 

0.379 

-80.469 

-0.697 

-0.821 

2 

80.409 

-0.849 

0.821 

-80.409 

-0.962 

-0.743 

3 

80.294 

-1.003 

0.743 

-80.294 

-0.769 

-1.039 

4 

80.123 

-0.792 

1.039 

-80.123 

-0.976 

-0.945 

5 

79.910 

-0.947 

0.945 

-79.910 

-0.825 

-1.190 

6 

79.660 

-0.729 

1.190 

-79.660 

-1.114 

-0.959 

7 

79.404 

-0.654 

0.959 

-79.404 

-1.267 

-0.529 

8 

79.178 

-0.376 

0.529 

-79.178 

-1.619 

0.516 

9 

79.004 

-0.188 

-0.516 

-79.004 

-1.757 

1.898 

10 

78.850 

0.690 

-1.898 

-78.850 

-2.349 

4.807 

11 

78.573 

0.666 

-4.807 

-78.573 

-1.367 

6.801 

12 

78.035 

0.058 

-6.801 

-78.035 

-0.058 

6.917 

13 

77.184 

-1.023 

-6.917 

-77.184 

1.023 

4.853 

14 

76.145 

-2.037 

-4.853 

-76.145 

2.037 

0.792 

15 

75.023 

-1.031 

-0.792 

-75.023 

1.031 

-1.300 

16 

74.050 

-1.584 

1.300 

-74.050 

1.584 

-4.481 

17 

73.321 

-1.250 

4.481 

-73.321 

1.250 

-7.004 

18 

72.923 

-0.676 

7.004 

-72.923 

0.676 

-8.361 

LC-5  CONSTRUCTION  CONDITION 
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Table  B5  (Continued) 


NUMBER  OF 

NODAL  POINT 

LOAD  CARDS  = 

19 

SPECIFIED  NOCAL  POINT 

LOADS 

NODE 

FX 

FY 

MOM 

1 

12.080 

0. 

2 

-4.25tJ 

23.770 

0. 

3 

-8.34k) 

22»6ii)id 

0. 

4 

-12.140 

20.710 

0. 

5 

-15.500 

18.160 

0. 

6 

-18.320 

15.060 

0. 

7 

-20.520 

11.520 

0. 

8 

-22.030 

7.700 

0. 

-22.830 

3.700 

0. 

Iw 

-22.900 

-0.360 

0. 

11 

-22.280 

-4.460 

0. 

12 

-21.000 

-8.350 

0. 

13 

-19.140 

-11.950 

0. 

14 

-16.760 

-15.130 

0. 

15 

-13.940 

-17.820 

0. 

16 

-10.760 

-19.960 

0. 

17 

-7.760 

-21.520 

0. 

18 

-3.700 

-22.470 

0. 

ly 

0. 

-11.470 

0. 

EliiMEWi  OJALb 

liJAD  TXPE  *  =1  CONCEM'RAIED  LOAD 

^2  UNIFORM  DISIHIBUI'ED  LOAD 
=3  CONCLNiHATED  MOMENT 
=4  TRIANGUIAR  DISTRIBUTED  LOAD 

LOCATION  CODE-  *LCAD  TYPES  1  AND  3: 

DISTANCE  FROM  I  NODE 
♦LOAD  TYPE  2;  =4) 

♦LOAD  TYPE  4: 

»ti  NODE  I  WEIGHTED,  -1  NODE  J  WEIGHTED 

DIRECTION  CODE-  ♦LOAD  TYPES  1-4: 

X=2,  Y=1 

ElEMENI  NODE  LOAD  MAGNITUDE  liXATION  DIRECTION 
I  J  TYPE  OF  LOAD  CODE  CODE 


(Continued) 
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Table  B5  (Continued) 


ElEMEMS  ON  AN  ELASTIC  FOUNCATTON 
EJJSMENl'  SPRING  PRESSURE 


1  0. 

2  0. 

3  0. 

4  0. 

5  0. 

6  0. 

7  0. 

8  0. 

9  0. 

10  0.502813E 

11  0. 82094 2E 

12  0.860922E 

13  0.732114E 

14  0. 52071 0E 

15  0.267984E 

16  0.130/8BE 

17  0. 

18  0. 


NUMBER  OF  ri'ERATTONS 


0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

00  0.104907E  00 

00  0.255364E  00 

00  0.224433E  00 

00  0.100573E  00 

00  0. 

00  U. 

00  0. 

0. 

0. 


O 


(Continued) 
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Table  B5  (Continued) 


V 


ELEMENl  ZCM, 

UNITS 

OF  INTEGRATION  (MAX*20) 

1  EVEN 

0. 

9 

2.008, 

2  EVEN 

9 

2.018, 

3  EVEN 

kJ. 

9 

2.009, 

4  EVEN 

9 

2.011, 

5  EVEN 

id. 

9 

2.008, 

6  EVEN 

ki. 

9 

2.017, 

7  EVEN 

d. 

9 

2.009, 

8  EVEN 

d. 

9 

2.008, 

9  EVEN 

k). 

9 

2.017, 

10  UDL> 

fc). 

9 

2.017, 

11  OD£> 

0. 

t 

2.008, 

12  ODD 

0. 

9 

2.009, 

13  ODD 

0. 

9 

2.017, 

14  ODD 

0. 

9 

1.909, 

1.994, 

15  ODD 

0. 

9 

1.341, 

2.029, 

16  ODD 

id. 

9 

0.836, 

2.009, 

17  EVEN 

d. 

9 

2.018, 

18  EVEN 

id. 

9 

2.007, 

NOCftL  tOlNT'  DISPLACfcMEWTS 
NODE  X-DISPL  Y-DISPL  ROT-DISPL 


1  i). 

2  -0.47287E-M3 

3  -i().818k)4E-B3 

4  -i3.93614E-83 

5  -0.8B432E-83 

6  -W.47999E-03 

7  -0.85794E-04 

8  0.21656E-03 

9  ifl.29724E-«)3 

10  0.11745E-03 

11  -0.23800E-03 

12  -0.65940E-03 

13  -0.10332E-02 

14  -0.12746E-02 

15  -0.13193E-02 

16  -0.11561E-02 

17  -0.85687E-03 

18  -0.45762E-03 

19  0. 


0. 10936E-01 
0.10724E-01 
0.10108E-01 
0. 92025E-02 
0.81554E-02 
0.71385E-02 
0.62605E-02 
0.55797E-02 
0.50529E-02 
0.45776E-02 
0.40588E-02 
0.34464E-02 
0.27425E-02 
0.19879E-02 
0.12730E-02 
ttf.68427E-03 
0.21872E-03 
-0.10786E-03 
-0.23032E-03 


0. 

-0.16994E-03 
-0.32487E'03 
-0. 43463E>03 
Hi).48074EH93 
-0.45301E-03 
-0.35482E-03 
-0.19755E-03 
-0.14S28E-04 
0.14198E-03 
0.24741E-03 
0. 30700E-03 
0.32063EHa3 
0.30040B-03 
0.22964EH33 
0.16238E-03 
0.12923E-03 
0.77802E-04 
0. 
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Table  B5  (Concluded) 


END  ACTIONS 

WITH  NODAL  POINT 

IDAC6 

eiEMEMl' 

AXIAL  I 

SHEAR  I 

HCMENT  I 

AXIAL  J 

SHEAR  J 

MQMENl'  J 

1 

136.791 

-0.182 

4.650 

-136.791 

0.182 

-5.016 

2 

136.786 

0.626 

5.016 

-136.786 

-0.626 

-3.751 

3 

136.782 

0.629 

3.751 

-136.782 

-0.629 

-2.489 

4 

136. 759 

1.173 

2.489 

-136.759 

-1.173 

-0.130 

5 

136.716 

0.915 

0.130 

-136.716 

-0.915 

1.707 

6 

136.625 

1.064 

-1.707 

-136.625 

-1.064 

3.853 

7 

136.459 

0.615 

-3.853 

-136.459 

-0.615 

5.089 

8 

136.215 

0.114 

-5.089 

-136.215 

-0.114 

5.318 

9 

135.857 

-0.880 

-5.318 

-135.857 

0.880 

3.543 

Id 

135.351 

-0.836 

-3.543 

-135.351 

0.187 

2.646 

11 

134.652 

-1.409 

-2.646 

-134.652 

0.308 

1.112 

12 

133.861 

-1.193 

-1.112 

-133.861 

0.098 

0.030 

13 

133.012 

-0.982 

-0.030 

-133.012 

0.149 

-0.896 

14 

132.155 

-1.384 

0.896 

-132.155 

0.907 

-3.000 

15 

131.350 

0.982 

3.000 

-131.350 

-1.153 

-0.736 

16 

130.649 

-0.247 

0.736 

-130.649 

0.193 

-1.138 

17 

130.543 

-0.315 

1.138 

-130.543 

0.315 

-1.772 

18 

130.259 

-0.439 

1.772 

-130.259 

0.439 

-2.654 

***  $END  *** 
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V 


I 


W  =  ROCK  LOAD,  tsf 

R.  =  INTERNAL  RADIUS 

R„  =  EXTERNAL  RADIUS 
e 

e  =  ANGLE  FROM  VERTICAL 

^  =  ANGLE  OF  SEGMENT 
X. 


=  HORIZONTAL  PROJECTION 
®  OF  EXTERNAL  FACE 

F  =  VERTICAL  ROCK  FORCE 


K  "  sin  (e  +  f)  -  R  sin  (e  -  |) 


Figure  B5.  Formula  for  calculating  rock  loads  for  example  2 


B37 


LOAD  CASE  1  =  ROCK  LOAD  +  WEIGHT  OF  LINING 

Rock  Load  =  0.275  tons/ft 

For  a  13-in.  lining,  =  12.083  ft 

0  =  30  degrees 

(j)  =  10  degrees 


At  Node  16, 

X  =  R  sin  (0  +  I-)  -  R  sin  (0  -  4) 
e  e  2'  e  '2' 

=  1.824  ft 

So,  if  F  =  WXg  ,  F  =  1.0032  kips  due  to  rock  load 
To  this,  add  the  load  due  to  the  weight  of  the  tunnel  lining: 
Weight  of  a  13-in.  lining  =  (tvR^  -  irR^)  0.150  kef 

=  11.78  kips/ft 

Therefore,  for  one  10-degree  segment,  the  weight  of  the  tunnel 
lining 

11.78  kips/ft 
36  segments 

=  0.32  kips/ft 


Total  Load  at  Node  16 

Add  F  to  the  weight  of  the  tunnel  lining: 

1.0032  kips  +  0.32  kips  =  1.32  kips 


Figure  B6.  Sample  calculation  of  nodal  loads  at  node  16  for  load  case  1, 

example  2 
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=  INTERNAL  WATER  HEAD  AT  INVERT,  ft 
^  EXTERNAL  WATER  HEAD  AT  INVERT,  ft 
=  INTERNAL  RADIUS 
=  EXTERNAL  RADIUS 
=  ANGLE  FROM  VERTICAL 
=  ANGLE  OF  SEGMENT 


^  FORCE  DUE  TO  INTERAL  HEAD 

=  FORCE  DUE  TO  INTERNAL  HEAD 
=  FORCE  DUE  TO  EXTERNAL  HEAD 
=  INTERNAL  SEGMENT  LENGTH 
=  EXTERNAL  SEGMENT  LENGTH 
=  NET  FORCE 


=  [(h.  -  R.)  .  cos  0R.]  S.Y^ 

F  =  [(h  -  R.)  -  cos  0R  ]  S  Y 

e  **  e  V  e  e*w 

irZR. 

c  =  - ] - 

i  No.  of  segments 

if2R 

5  s  - § - 

e  No.  of  segments 


Figure  B7.  Formula  for  calculating  water  loads  for  example  2 


LOAD  CASE  2  =  ROCK  LOAD  +  WEIGHT  OF  LINING  +  INTERNAL  HEAD 


For  a  13-in.  lining,  36-segment  model,  ^  =  10  degrees 
R.  =  11  ft, 

so  S.  =  =  1.920  ft 

Rg  =  12,083  ft, 

so  S  =  ?li-12.083)  ^  2.109  ft 

e  36 

For  Node  16, 

F.  =  [(h.  -  R.)  -  R.  cos  0]  S.Y^^ 

=  [(180  -  11)  -  11  cos  30](1. 92)0. 0624 
=  19.11  kips 

F^  =  -  R,-)  -  R„  cos  e]  S  Y 

e  e  1  e  e  w 

=  [(150  -  11)  -  12.083  cos  30](2. 109)0. 0624 

=  16.92  kips 

=  16.92  -  19.11  =  -2.19  kips 
F^  =  F^  sin  0  =  -2.19  sin  30  =  1.10  kips-*- 

Fy  =  F^  cos  0  =  -2.19  cos  30  =  1.897  kipst 

Adjust  Fy 

The  Rock  Load  and  Lining  Weight  from  Load  Case  1  was  1.32  kips. 
So, 

Fy  =  1.897  -  1.32  =  0.58  kipst 


Figure  B8.  Sample  calculation  of  nodel  loads  at  node  16  for  load  case  2, 

example  2 


LOAD  CASE  3  =  ROCK  LOAD  +  WEIGHT  OF  LINING  +  EXTERNAL  HEAD 
For  Node  16, 

F.  =  [(100  -  11)  -  11  cos  30](1. 92)0. 0624 
=  9.52  kips 

Fg  =  [(185  -  11)  -  12.083  cos  30](2. 109)0.0624 

=  21.52  kips 

F  =  F  -  F. 
n  e  1 

=  21.52  -  9.52  =  12.00  kips 
F^  =  12.0  sin  30  =  -6.0  kips-t- 
Fy  =  12.0  cos  30  =  -10.39  kips+ 

Adjust  F 

The  Rock  Load  and  Lining  Weight  from  Load  Case  1  was  1.32  kips 
So, 

Fy  =  -10.39  +  (-1.32)  =  -11.71  kipsf 


Figure  B9.  Sample  calculation  of  nodal  loads  at  node  16  for  load  case  3, 

example  2 


JOINT  LOCATION 

b.  Load  case  2 

moment  in  the  concrete  lining  for  example 


JOINT  LOCATION 

Load  case  5 


In  accordance  with  letter  from  DAla-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Pacsiisile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a  facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 


Holtham,  William  J. 

User’s  guide  :  Computer  program  for  design  and  analysis 
of  cast-in-place  tunnel  linings  (NEWTUN)  :  final  report  / 
by  William  J.  Holtham,  James  R.  Fay,  U.S,  Army  Engineer 
Division,  New  England  ;  prepared  for  Office,  Chief  of 
Engineers,  U.S,  Army  ;  monitored  by  Automatic  Data 
Processing  Center,  U.S.  An^y  Engineer  Waterways  Experiment 
Station.  —  Vicksburg,  Miss.  ;  U.S.  Army  Engineer  Waterways 
Experiment  Station  ;  Springfield,  Va.  :  available  from 
NTIS,  1981. 

15,  [i*6]  p.  :  ill.  ;  27  cm.  —  (Instruction  report  /U.S. 
An^y  Engineer  Watei^mys  Experiment  Station  ;  K-8I-U) 

Cover  title. 

"March  1981." 

"A  report  under  the  Computer-Aided  Structural  Engineering 
(CASE)  Project." 

Bibliography:  p.  15. 

1.  Mathematical  analysis.  2.  NEWTUN  (Computer  program). 
3-  Tunnels.  I.  Pay,  James  R.  II,  United  States. 

Army.  Corps  of  Engineers.  New  England.  III.  United 
States.  Army.  Corps  of  Engineers.  Office  of  the  Chief  of 


Holtham,  William  J. 

User’s  guide  :  Computer  program  for  design  :  ...  I98I. 

(Card  2) 

Engineers.  IV.  United  States.  Army  Engineer  Waterways 
Experiment  Station.  Automatic  Data  Processing  Center. 

V.  Title  VI.  Series:  Instruction  report  (United 

States.  Army  Engineer  Waterways  Experiment  Station)  ;  K-81-U. 

TA7.W3UI  no.K-81-U 


WATERWAYS  EXPERIMENT  STATION  REPORTS 
PUBLISHED  UNDER  THE  COMPUTER-AIDED 
STRUCTURAL  ENGINEERING  (CASE)  PROJECT 


Technical  Report  K-78-1 

Instruction  Report  0-79-2 

Technical  Report  K-80-1 
Technical  Report  K-80-2 

Instruction  Report  K-80-1 

Instruction  Report  K-80-3 

Instruction  Report  K-80-4 

Instruchon  Report  K-80-6 

Instruction  Report  K-80-7 

Technical  Report  K-80-4 

Technical  Report  K-80-5 
Instruction  Report  K-81-2 

Instruction  Report  K-81-3 

instruction  Report  K-81-4 
Instruction  Report  K-81-6 

Instruction  Report  K-81-7 


Title 

List  of  Computer  Programs  for  Computer-Aided 
Structural  Engineering 

User's  Guide:  Computer  Program  with  Interactive 
Graphics  for  Analysis  of  Plane  Frame  Structures 
(CFRAME) 

Survey  of  Bridge-Oriented  Design  Software 

Evaluation  of  Computer  Programs  for  the 
Design/Analysis  of  Highway  and  Railway  Bridges 

User's  Guide;  Computer  Program  for  Design/Review 
of  Curvilinear  Conduits/Culverts  (CURCON) 

A  Three-Dimensional  Finite  Element  Data 
Edit  Program 

A  Three-Dimensional  Stability  Analysis^Design 
Program  (3DSAD) 

Report  1:  General  Geometry  Module 

Basic  User’s  Guide  Computer  Program  for  Design  and 
Analysis  of  Inverted-T  Retaining  Walls  and  Floodwalls 
(TWDA) 

User’s  Reference  Manual;  Computer  Program  for 
Design  and  Analysis  of  Inverted-T  Retaining  Walls 
and  Floodwalls  (TWDA) 

Documentation  of  Finite  Element  Analyses 
Report  1;  Longview  Outlet  Works  Conduit 
Report  2;  Anchored  Wall  Monolith.  Bay 
Springs  Lock 

Basic  Pile  Group  Behavior 

User's  Guide:  Computer  Program  for  Design  and 
Analysis  of  Sheet  Pile  Walls  by  Classical  Methods 
(CSHTWAL) 

Report  1;  Computational  Processes 
Report  2;  Interactive  Graphics  Options 

Validation  Report:  Computer  Program  for  Design 
and  Analysis  of  Inverted-T  Retaining  Wails  and 
Floodwalls  (TWDA) 

User's  Guide:  Compu  er  Program  for  Design  and 
Analysis  of  Cast-in-Place  Tunnel  Linings  (NEWTUN) 

User's  Guide:  Computer  Program  for  Optimum 
Nonlinear  Dynamic  Design  of  Reinforced  Concrete 
Slabs  Under  Blast  Loading  (CBARCS) 

User's  Guide:  Computer  Program  for  Design  or 
Investigation  of  Orthogonal  Culverts  (CORTCUL) 


Date 
Feb  1978 

Mar  1979 

Jan  1980 
Jan  1980 

Feb  1980 

Mar  1980 

Jun  1980 
Dec  1980 

Dec  1980 

Dec  1980 
Dec  1980 

Dec  1980 

Feb  1981 
Mar  1981 

Feb  1981 

Mar  1981 

Mar  1981 

Mar  1981 


